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Executive  Summary 


The  response  of  a  supersonic  high  Reynolds  number  turbulent  boundary  layer  flow  subjected 
to  mechanical  distortions  was  experimentally  examined.  Local  disturbances  were  introduced  via 
small-scale  wall  patterns,  and  global  distortions  were  induced  through  streamline  curvature 
driven  pressure  gradients.  Local  surface  topologies  included  A:-type  diamond  and  <i-type  square 
elements;  a  smooth  wall  was  examined  for  comparison  purposes.  Three  global  distortions  were 
studied  with  each  of  the  three  surface  topologies.  Measurements  included  planar  contours  of  the 
mean  and  fluctuating  velocity  pressure  sensitive  paint  and  schlieren  photography.  The  velocity 
data  were  acquired  with  sufficient  resolution  to  characterize  the  mean  and  turbulent  flow 
structure  and  to  examine  interactions  between  the  local  surface  roughness  distortions  and  the 
imposed  pressure  gradients  on  the  turbulence  production.  A  strong  response  to  both  the  local  and 
global  distortions  was  observed  with  the  diamond  elements,  where  the  effect  of  the  elements 
extended  into  the  outer-regions  of  the  boundary  layer.  It  was  shown  that  the  primary  cause  for 
the  observed  response  was  the  result  of  local  shock  and  expansion  waves  modifying  the 
turbulence  structure  and  production.  By  contrast,  the  square  elements  showed  a  less  pronounced 
response  to  local  flow  distortions  as  the  waves  were  significantly  weaker.  However,  the  frictional 
losses  were  higher  for  the  blunter  square  roughness  elements.  Detailed  quantitative 
characterizations  of  the  turbulence  flow  structure  and  the  associated  production  mechanisms  are 
described  herein.  These  experiments  demonstrate  fundamental  differences  between  supersonic 
and  subsonic  rough-wall  flows,  and  the  new  understanding  of  the  underlying  mechanisms 
provides  a  scientific  basis  to  systematically  modify  the  mean  and  turbulence  flow  structure  all 
the  way  across  supersonic  boundary  layers. 

This  report  was  based  on  the  article  “Response  of  supersonic  turbulent  boundary  layers  to 
local  and  global  mechanical  distortions,”  by  Ekoto  et  al  in  the  Journal  of  Fluid  Mechanics,  vol. 
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630,  2009,  pp.  225-265. 


Nomenclature 


Ab  =  surface  element  blockage  area 

a  =  roughness  element  width 

B,  =  [(T„-rj/Tj'- 

B2  = 

B;  =  [4B‘+B;f- 

b  =  roughness  length 

d  =  distortion  parameter  [  {du Idx  +  dvldx  +  dvj dy)  /  {du  jdy)  ] 

k  =  roughness  height 

=  turbulent  kinetic  energy  {-‘A tI  ! p) 

M  =  Mach  number 

p  =  pressure 

Apw  =  axial  wall  pressure  variation  [{p„  —  p„)lp„] 

Pij  =  turbulent  stress  production  ( ^  ) 

P  =  turbulent  kinetic  energy  production 

r  =  turbulent  recovery  factor  (=0.89) 

Re  =  Reynolds  number  (puL/p) 

Sxy  =  Strain  Rate  [=  j  (du  /dy  +  dv  /  dx)  ] 

T  =  temperature 

Taw  =  adiabatic  wall  temperature  (=T^{\  +  \r(Y-  1)M^ ] ) 
u,v  =  axial  and  transverse  velocity  components 
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=  van  Driest  II  scaled  velocity  (=  [sin  '  { [25j^  (w  /  }  +  sin  '  {B^  I B^)]  I B^) 

=  reference  nozzle  exit  velocity  (=  625  m/s) 

1  /9 

=  friction  velocity  (Tm//?w) 

=  roughness  spacing 

=  Cartesian  coordinates  (axial,  wall  normal  and  lateral  directions,  respectively) 

=  turbulent  stress  turning  angle 
=  boundary  layer  thickness 
=  boundary  layer  displacement  thickness 
=  reference  nozzle  exit  boundary  layer  thickness  (=10  mm) 

=  flow  turning  angle  relative  to  the  test  section  coordinate  system 
=  von  Karman  constant  (=  0.41) 

=  roughness  element  spacing  wavelength 
=  dynamic  viscosity 
=  density 

=  momentum  thickness 
=  Reynolds  shear  stress  {=-pu[u'j  ] 

=  wall  shear  stress 
=  Vorticity  [=  \  (9v  I  dx-du  I  dy)  ] 


c  =  izTpuJpy 

subscripts 

aw  =  Adiabatic  Wall  Condition 


index 


i,j,k 

L  =  laminar 

0  =  total  condition 

ref  =  reference  condition 

5  =  sand  grain 

w  =  wall 

A  =  turbulence  stress  principal  axis 

oo  =  Freestream  condition 


superscripts 


T  =  turbulent 


0  =  Reynolds  mean  component 

0'  =  Reynolds  fluctuation 

+  =  inner  variable  scaled  component  [p^u* {^/ fJ^w\ 
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1.  Introduction 


1.1  Research  Objective  and  Approach 

Improved  understanding  of  the  underlying  physical  processes  associated  with  mechanically 
distorted  turbulent  supersonic  boundary  layers  with  surface  roughness  has  many  applications 
including  reentry  vehicle  thermal  protection  and  internal  aero-propulsive  flows.  Relevant  surface 
topologies  range  from  random  pitting,  spalting  and  contamination  deposits,  to  periodic  elements 
associated  with  tiled  thermal  protection  systems  and  ablative  thermal  protection  systems  with  an 
underlying  honeycomb  structure.  In  high  Reynolds  number  supersonic  turbulent  boundary  layer 
flows,  the  sonic  point,  the  point  above  which  the  flow  is  supersonic,  is  close  to  the  surface.  Thus, 
moderate  sized  roughness  elements  may  extend  into  the  supersonic  portion  of  the  boundary 
layer,  causing  corresponding  shock  and  expansion  waves  that  can  alter  the  mean  and  turbulence 
flow  structure  in  a  manner  qualitatively  different  than  in  subsonic  flows. 

The  objective  of  the  present  study  was  to  achieve  improved  understanding  of  the  underlying 
physical  processes  associated  with  mechanically  distorted  turbulent  supersonic  boundary  layers 
with  surface  roughness.  To  meet  this  objective,  a  series  of  experiments  were  performed  in  a 
supersonic  (M  =  2.87)  high  Reynolds  number  boundary  layer  (Ree  =  40,000).  To  provide  a 
complete  and  relevant  basis  for  analyses,  experiments  were  performed  with  the  same 
instrumentation  for  the  following  cases:  (1)  zero  pressure  gradient  boundary  layer  on  a  smooth 
wall,  (2)  zero  pressure  gradient  boundary  layer  with  square,  canonical  J-type,  roughness 
topology,  (3)  zero  pressure  gradient  boundary  layer  with  diamond  shaped  A:-type  roughness 
topology,  which  produces  a  strong  secondary  flow  structure  through  a  pattern  of  attached  oblique 
shocks  and  expansion  waves,  (4)  favorable  pressure  gradient  boundary  layers  on  a  smooth  wall, 
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(5)  favorable  pressure  gradient  boundary  layers  with  square  roughness,  (6)  favorable  pressure 
gradient  boundary  layers  with  diamond  roughness,  (7)  combined  favorable  then  adverse  pressure 
gradient  boundary  layer  on  a  smooth  wall,  (8)  combined  pressure  gradient  boundary  layer  with 
square  roughness,  and  (9)  combined  pressure  gradient  boundary  layer  with  diamond  roughness. 
The  first  two  cases  provide  a  canonical  background  for  comparison  purposes,  where  the  t/-type 
roughness  was  expected  to  follow  subsonic  trends.  Case  three  characterized  interactions  between 
the  local  pressure  waves  and  boundary  layer  mean  and  turbulence  flow  structure.  Cases  four,  five 
and  seven  through  nine  were  examined  to  characterize  global  distortion  effects  for  comparison  to 
the  local  gradients  in  case  3.  The  favorable  pressure  gradients  in  case  six  were  selected  isolate 
and  demonstrate  the  importance  of  shock-turbulence  interactions  near  the  leading  edge  of  the 
roughness  elements.  The  experiments  provided  quantification  of  the  effects  of  mechanical 
distortions  on  the  mean  flow  properties  and  turbulence  statistics,  structure  and  production. 

This  study  builds  on  three  basic  bodies  of  knowledge:  (1)  supersonic  rough  wall  boundary 
layers,  (2)  distorted  supersonic  turbulent  boundary  layers,  and  (3)  shock  wave  turbulence 
interactions.  The  review  presented  below  summarizes  the  key  building  blocks  for  the  present 
study. 

1.2  Supersonic  Rough  Wall  Boundary  Layers 

The  current  understanding  of  high-speed  rough  wall  mean  and  surface  flow  processes  are 
founded  to  a  large  degree  in  applying  compressibility  scaling  to  the  corresponding  low-speed 
incompressible  database  [e.g.,  see  Liepmann  and  Goddard  (1957),  Goddard  (1959),  Morkovin 
(1961),  Reda  (1975),  Berg  (1979),  Latin  and  Bowersox  (2000,  2002)].  Incompressible  rough 
wall  boundary  layers  have  been  the  subject  of  considerable  attention;  e.g.,  see  Nikuradse  (1933), 
Perry  et  al.  (1969,  1987),  Schlichting  (1979),  Ligrani  and  Moffat  (1986),  Jimenez  (2004), 
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Shockling  et  al.  (2006)  and  Schultz  and  Flack  (2007).  These  studies  have  shown  that  surfaee 
roughness  has  a  direct  effect  on  the  inner-region  of  the  law  of  the  wall  and  is  often  deseribed  by 
the  single  roughness  Reynolds  number  parameter.  Nikuradse  (1933)  demonstrated  that  for  low- 
speed  flows,  sand-grain  generated  roughness  increased  the  velocity  defect  and  skin  friction,  and 
shifted  the  veloeity  law  of  the  wall  plot  downward.  The  amount  of  downward  shift  was  shown  to 
be  a  function  of  the  roughness  height  Reynolds  number  k'^ .  Nikuradse’s  data  (1933)  showed  no 
measurable  effect  of  roughness  for  less  than  about  5-7,  a  condition  commonly  referred  to  as  a 
hydraulieally  smooth.  On  the  other  hand,  for  k^  greater  than  about  60  -  70,  the  frictional  losses 
and  the  velocity  profile  were  independent  of  viscosity,  and  the  flow  is  ealled  fully  rough.  In  the 
transition  regime  where  1  <  k^  <  70,  the  flow  properties  depend  on  both  the  viscosity  and  the 
roughness.  The  above  descriptions  were  developed  for  tightly  paeked  sand  grain  roughness. 
Sehliehting  (1979)  developed  an  equivalent  sand  grain  roughness  to  relate  other  roughness 
topologies  to  the  sand  grain  roughness.  The  relation  between  the  aetual  roughness  height  and  the 
equivalent  sand  grain  height  is  Ink^ /A:  =  /ir(8.5-C),  where  C  is  the  intereept  for  the  aetual 

roughness  (Aw^  =  InA:^ -1-5-C).  Generally,  the  location  of  the  profile  origin  is  selected 

somewhere  between  the  base  and  top  of  the  elements  to  improve  the  quality  of  the  veloeity 
profile  in  the  logarithmic  layer  [Perry  et  al.  (1969)]. 

Goddard  (1959)  demonstrated  that  the  sand  grain  roughness  velocity  shift  for  high-speed 
flow  eollapsed  onto  the  incompressible  eorrelation  when  van  Driest  II  (1951)  scaling  was 
employed;  that  is  Ink^  -1-5-C,  where  Me# is  the  veloeity  with  van  Driest  II  sealing. 

Goddard  (1959)  also  found  that,  for  fiilly  rough  sand  grain  flow,  C  matched  the  low-speed  value 
of  8.5  with  van  Driest  II  scaling.  The  eompilation  of  the  velocity  shift  data  shown  in  Fig.  la 
confirms  Goddard’s  conclusion. 
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The  equivalent  sand  grain  description  has  proven  inadequate  to  characterize  roughness 
topology  effects.  Hence,  investigations  of  the  influence  of  the  roughness  geometry  have  also 
been  performed  [Perry  et  a.  (1969),  Simpson  (1973),  Antonia  and  Wood  (1975),  Rons  et  al. 
(2001),  Jimenez  (2004),  McClain  et  al.  (2003)].  Overall,  a  general  description  has  not  been 
developed.  However,  two  categories  of  roughness  elements  have  been  defined:  k-type  and  d-type. 
For  k-type  roughness  elements,  the  equivalent  sand  grain  height  is  proportional  to  the  roughness 
height,  where  for  d-type,  the  equivalent  sand  grain  height  is  proportional  to  the  boundary  layer 
height  [Jimenez  (2004)].  The  solidity  [Schlichting  (1979)]  is  often  used  as  a  correlating 
parameter.  However,  different  definitions  have  been  proposed  in  the  articles  listed  above. 

Raupach  et  al.  (1991)  provides  a  comprehensive  survey  (157  articles)  of  low-speed  scaling 
for  laboratory  and  atmospheric  scale  flow  fields.  They  extended  the  Reynolds  number  similarity 
concept  of  Townsend  (1976)  to  include  rough  walls  by  defining  wall  similarity  as  the  boundary- 
layer  turbulent  motions  outside  the  roughness  or  viscous  sublayer  at  high  Reynolds  numbers  that 
are  independent  of  wall  roughness  and  viscosity.  The  exceptions  were  the  role  of  the  wall  in 
setting  the  friction  velocity,  effective  origin  and  boundary-layer  thickness.  Raupach  et  al.  also 
pointed  out  that  an  important  consequence  of  the  wall  similarity  hypothesis  is  that  the  organized 
structures  (double-roller  eddy  or  A- vortex)  in  the  outer-region  of  the  rough  wall  boundary  layers 
should  be  the  same  as  those  for  smooth  wall  boundary  layers,  even  though  the  near  wall  vortex 
generation  for  smooth  and  rough  walls  is  fundamentally  different.  Smooth  walls  are 
characterized  by  the  bursting  process  and  rough  surfaces  by  organized  wake  vortices.  For  most 
roughness  types,  the  near  wall  flow  is  so  disturbed  that  the  smooth  wall  low-speed  streaks  are 
eradicated;  however,  there  is  evidence  of  a  rough  wall  counterpart  [Raupach  et  al.  (1991)].  Grass 
(1971)  performed  detailed  investigations  of  the  influence  of  surface  roughness  on  the  turbulent 
structure  and  turbulence  production.  His  data  confirmed  that  the  smooth  wall  inrush  and  ejection 
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process,  typically  associated  with  the  bursting  process,  correlate  with  a  high  contribution  to  the 
Reynolds  stress  and  hence  turbulence  production  close  to  the  wall.  The  inrush  and  ejection 
process  was  found  to  be  present  irrespective  of  the  surface  roughness  condition.  Grass  also  noted 
a  strong  interaction  between  the  inner  and  outer  flows.  In  a  recent  experimental  study,  George 
and  Simpson  (2000)  resolved  the  flow  around  the  roughness  elements  and  obtained  data  very 
close  to  the  wall.  These  data  show  that  the  near  wall  structure  was  highly  dependent  on  the 
roughness  topology.  However,  the  outer-region,  when  scaled  by  the  friction  velocity,  was 
independent  of  the  roughness,  which  further  confirmed  the  conclusions  of  Raupach  et  al.  (1991). 
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Fig.  1  High-speed  mean  flow  properties  [Theory:  B  =  5.5,  C  =  8.5,  K=  0.4] 


High-Speed  studies  that  include  turbulence  data  are  less  common.  The  high-speed  rough  wall 
Reynolds  shear  stress  measurements  of  Latin  and  Bowersox  (2000)  are  compared  to 
representative  low-speed  data  along  with  additional  high-speed  smooth  wall  data  in  Fig.  lb.  In 

Fig.  lb,  =-puv  and  =  pJJ^  .  The  laminar  term,  is  only  significant  near  the  wall 
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and  in  the  logarithmic  region  of  the  boundary  layer,  The  scatter  in  the  data  in  this  plot  is 

representative  of  the  current  database  [e.g.,  see  the  compilations  in  Dussauge  et  al.  (1996)].  The 
supersonic  smooth  wall  data  of  Luker  et  al.  (2000)  and  the  supersonic  sand-grain  and  square 
roughness  topology  roughness  data  of  Latin  and  Bowersox  (2000)  also  followed  Townsend’s 
(1976)  scaling,  with  Morkovin’s  p  /  p„  scaling  (1961).  However,  the  2-D  bar  roughness  element 

topology  studied  by  Latin  and  Bowersox  (2000)  did  not  follow  the  pattern.  These  data  indicated 
a  dependence  of  outer-region  turbulent  levels  on  the  roughness  topology.  Latin  and  Bowersox 
(2000)  also  found  that  the  correlation  between  transverse  velocity  and  density  fluctuations  ( 

p'v  /  pu ),  directly  measured  by  combining  the  response  from  hot-wire  anemometry  with  laser 
Doppler  velocimetry,  demonstrated  a  systematic  dependence  on  the  roughness  height.  This  was 
the  general  conclusion  for  all  of  the  turbulence  properties  that  included  a  thermodynamic 
variable.  They  believed  this  dependence  was  the  result  of  the  roughness  elements  generated 
shock  and  expansion  waves  interacting  with  the  boundary  layer  turbulence.  The  present  study 
showed  that  secondary  distortions  associated  with  such  waves  significantly  affect  the  transport  of 
turbulence. 

1.3  Distorted  Supersonic  Turbulent  Boundary  Layers 

Curvature  driven  bulk  compression  and  dilatation  is  important  to  the  present  study  on  both 
the  local  and  global  scale.  Considerable  theoretical  and  experimental  studies  have  been 
performed  to  characterize  the  response  of  supersonic  boundary  layers  to  bulk 
compression/dilatation  [e.g.,  see  Bradshaw  (1969,  1974),  Dussauge  and  Gaviglio  (1987),  Smits 
et  al.  (1989),  Smith  and  Smits  (1991),  Smits  and  Dussauge  (1996),  Spina  et  al.  (1994),  Amette  et 
al.  (2000)  and  Luker  et  al.  (2000)].  The  data  compiled  in  Spina  et  al.  (1994)  show  that  adverse 
pressure  gradients  and/or  concave  wall  curvature  result  in  destabilized  flow  with  amplified 
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turbulence  levels,  where  favorable  pressure-gradients  and  convex  wall  curvature  have  a 
stabilizing  effect.  Bradshaw’s  (1974)  distortion  {d)  parameter,  defined  as  the  ratio  of  the 
secondary  strain-rates  to  the  primary  velocity  gradient  normal  to  the  wall,  is  often  used  to 
characterize  the  strength  of  the  pressure  gradient.  A  distortion  is  generally  considered  mild  if  <7  < 
0.01  and  strong  for  <7  >  0.1.  Bradshaw  also  pointed  out  that  the  effect  of  the  extra  strain-rates  on 
the  turbulence  properties  is  an  order  of  magnitude  stronger  than  implied  by  the  extra  production 
alone.  Spina  et  al.  (1994)  utilized  two  impulse  parameters  to  characterize  the  imposed  distortion 
strength:  the  bulk  dilatation  impulse  parameter,  defined  as  In  p^l  p^,  where  subscripts 

“1”  and  “2”  correspond  to  the  pressure  before  and  after  the  imposed  pressure  gradient,  and  the 
streamline  curvature  parameter  given  by  =  A^,  where  (p  is  the  change  in  the  wall  angle  in 

radians.  Imposed  impulses  of  Ip  =  0.5  and  0.8  produced  200%  and  600%  increases  in  the  axial 
turbulent  shear  stress  intensities,  respectively.  Wall  curvature  was  found  to  further  enhanced  the 
amplifications.  The  studies  of  Amette  et  al.  (2000)  and  Luker  et  al.  (2000)  demonstrated  that 
axial  turbulence  intensity  reductions  of  up  to  90%  were  possible  for  imposed  wall  induced 
favorable  pressure  gradients.  Near  the  wall,  turbulent  shear  stress  reductions  of  up  to  75%  were 
observed,  and  in  the  outer  half  of  the  boundary  layer,  the  turbulent  shear  stresses  were  negative, 
while  the  principal  strain-rate  remained  positive.  Luker  et  al.  (2000)  experimentally 
demonstrated  that  the  associated  extra  production  terms  due  to  the  induced  strain-rates  resulted 
in  overall  negative  turbulence  production  in  the  same  region  of  the  boundary  layer.  Luker  et  al. 
(2000)  also  inferred  from  their  data  that  the  favorable  pressure  gradient  flow  led  to  large-scale 
eddy  disintegration,  resulting  in  increased  turbulent  energy  available  for  flow  dissipation  and  the 
observed  stabilizing  effect  on  the  boundary  layer.  The  present  study  demonstrated  that  the  trends 
cited  for  the  global  distortions  also  hold  for  small-scale  local  distortions  in  supersonic  flows. 
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1.4  Shock  Turbulence  Interaction 


Shock  turbulence  interaetion  is  important  to  the  present  study  on  a  loeal  scale  as  the  present 
surfaee  topology  produeed  loeal  shoek  waves.  As  deseribed  in  Mahesh  et  al.  (1995),  shoek 
turbulenee  interaetion  has  numerous  practical  applications.  Hence,  considerable  theoretical, 
experimental  and  numerieal  researeh  has  been  performed  to  eharacterize  the  response  of 
turbulenee  as  it  passes  though  a  shoek  wave  [Keller  and  Merzkirch  (1990),  Dolling  (1993), 
Jaequin  et  al.  (1993),  Mahesh  et  al.  (1993,  1995,  1997),  Thivet  et  al.  (2001),  Sinha  et  al.  (2003, 
2005)].  One  of  the  major  interests  for  the  present  study  is  the  shear  stress;  i.e.,  the  anisotropie 
portion  of  the  stress  tensor.  Recent  computational  models  [Wallin  and  Johansson  (2000)]  relate 
the  standard  anisotropie  portion  of  the  Reynolds  stress  tensor  to  the  mean  strain-rate  tensor 
relation  augmented  by  higher  order  strain  and  rotation-rate  tensors.  Additionally,  Mahesh  et  al. 
(1993)  demonstrated  that  the  evolution  of  the  Reynolds  shear  stress  in  an  anisotropic  shear  layer 
was  strongly  dependent  on  the  angle  between  the  compression  and  shear  layer.  Sinee  the  stress 
tensor  ean  be  diagonalized  through  the  choice  of  an  appropriate  eoordinate  system,  i.e.,  the 
prineipal  axes,  there  is  no  inherent  distinction  between  the  normal  and  shear  stresses,  only 
between  the  isotropie  and  anisotropic  portions  of  the  stress  tensor  [Miles  (2007)],  and  the 
prineipal  axes  should  be  heavily  biased  towards  the  orientation  of  a  shock  wave.  The  important 
parameters  needed  to  deseribe  the  turbulence  stress  evolution  are  the  angle  {a)  necessary  to 
rotate  the  normal  stresses  to  their  principal  axes  and  the  turbulent  kinetic  energy  (k^),  whieh 
remains  constant  regardless  of  orientation.  Mahesh  et  al.  (1997)  demonstrated  that  the  sign  of  the 

eorrelation  between  the  axial  veloeity  fluetuation  and  the  temperature  fluctuations,  that  is  uT' 

had  an  important  impact  on  the  role  of  the  shoek  turbulenee  interaction.  If  uT'  <  0,  then  the 
amplification  of  the  turbulenee  levels  (kinematic  and  thermodynamic)  across  the  shock  is 
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enhanced.  This  is  the  case  for  adiabatic  supersonic  turbulent  boundary  layers.  However,  if  u'T'  > 
0  upstream  of  the  shock,  which  would  be  the  case  for  a  highly  cooled  wall,  then  the  amplification 
of  the  kinetic  energy  is  strongly  suppressed.  Mahesh  et  al.  (1997)  also  demonstrated  that  the 
phase  difference  between  the  upstream  vorticity  and  entropy  waves  determines  whether 
amplification  occurs.  It  was  demonstrated  in  the  present  study  that  roughness  element  shock 
waves  enhanced  the  turbulence  levels  through  this  interaction  mechanism. 

2.  Facilities  and  Diagnostic  Test  Equipment 


2.1  Wind  Tunnel 

The  experiments  were  performed  in  a  supersonic  wind  tunnel  located  at  Texas  A&M 
University.  A  schematic  of  the  wind  tunnel  is  given  in  Fig.  2.  The  facility  was  operated  in  a 
blow-down  mode,  with  a  maximum  nominal  run  time  of  20  minutes.  The  freestream  Mach 
number  at  the  nozzle  exit  was  2.87  with  a  ±1.0%  variation  across  nozzle  exit  core  flow.  The 
nozzle  exit  dimensions  were  7.62  cm  by  7.62  cm.  The  stagnation  pressure  and  temperature  were 
set  to  690  kPa  (±5.0%)  and  315K  (±3.5%),  respectively.  At  the  operating  stagnation  temperature, 
the  adiabatic  wall  temperature  was  approximately  equal  to  the  ambient  temperature  in  the 
laboratory,  which  when  coupled  with  the  relatively  low  thermal  conductivity  of  the  Plexiglas 
floors  in  the  tunnel  resulted  in  minimal  wall  heat  transfer. 

A  pitot  probe  with  an  Omega  brand  PX303  pressure  transducer  (0  -  690  kPa  and  ±0.25% 
full-scale  accuracy)  was  used  to  monitor  the  stagnation  pressure  within  the  settling  chamber.  The 
total  temperature  was  measured  using  an  Omega  brand  JQSS  thermocouple  amplified  by  an 
OMNI  AMP-IV  thermocouple  amplifier.  Tunnel  static  pressure  was  measured  with  five  Omega 
brand  PX303  pressure  transducers  (0  -  103  kPa  and  ±0.25%  full-scale  accuracy).  Data  from  the 
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pressure  transducers  and  thermocouples  were  collected  with  a  Keithley  brand  KPCI-3110  data 
acquisition  board  with  a  sampling  rate  of  1.25  MS/s.  The  board  had  12-bit  resolution. 

Measurements  took  place  in  the  floor  boundary  layer  of  the  wind  tunnel.  The  test  section  was 
located  just  downstream  of  the  nozzle  exit  (see  Fig.  2).  Due  to  the  half  block  design,  the  floor 
boundary  layer  experienced  20.5  cm  of  constant  pressure  flow  in  the  axial  direction  at  the 
freestream  conditions  before  reaching  the  test  section  entrance.  The  boundary  layer  thickness  at 
the  nozzle  exit  was  1.0  cm. 


Fig.  2  Schematic  of  supersonic  wind  tunnel  (flow  is  up  and  to  the  right). 

2.2  Roughness  topologies 

Two  roughness  topologies  were  designed  to  produce  fundamentally  different  near  wall 
distortions.  Smooth  wall  models  were  also  tested  for  comparison  purposes.  The  models  were 
manufactured  from  Plexiglas  and  polished  using  Mothers  05100  Mag  and  Aluminum  Polish.  The 
estimated  roughness  height,  k*,  for  the  smooth  models  was  0.14. 
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A  square  roughness  topology  was  designed  to  produee  the  eanonical  “sand-grain”  behavior 
observed  in  low-speed  flow.  More  speeifieally,  Latin  and  Bowersox  (2000)  demonstrated  that  the 
3-D  machined  square  topologies  followed  the  accepted  correlations  (see  Fig.  lb).  To  further 
foster  the  accepted  scaling,  the  roughness  elements  spacing  was  selected  to  produce  the  t/-type 
roughness  characteristic.  Perry  et  al.  (1969)  remarked  that  for  t/-type  periodic  roughness 
elements,  recirculation  is  set  up  in  the  grooves  and  eddy  shedding  from  the  elements  into  the 
flow  is  negligible.  For  the  zero  pressure  gradient  smooth  wall  model,  the  law-of-the-wall  with 
van  Driest  II  (1951)  scaling,  indicated  that  the  sonic  line  was  nominally  0.12  mm  above  the  wall. 
For  the  rough- wall,  it  was  desired  to  have  a  fully  rough  topology,  while  maintaining  the 
roughness  height  below  the  sonic  line  to  minimize  roughness  element  induced  shock-wave 
generation.  With  the  correlation  in  Fig.  la,  the  roughness  height  was  selected  at  0.8  mm.  At  this 
y-location,  the  local  Mach  number  was  estimated  at  0.9.  Thus,  strong  localized  shock  waves  were 
not  expected  for  this  roughness  topology.  Instead,  the  roughness  was  expected  to  produce 
enhanced  momentum  loss  similar  to  that  for  low-speed  flows.  The  model  was  machined  from 
Plexiglas.  The  channels  were  created  using  a  1.6  mm  diameter  ball  end  mill  bit  and  were  0.8  mm 
deep,  creating  roughness  elements  with  a  length  and  width  of  1.6  mm,  and  separated  from 
adjacent  elements  by  a  1.6  mm  wide  channel.  A  schematic  of  the  element  shape  is  given  in  Fig. 
3a,  while  the  expected  flow  structure  is  illustrated  in  Fig.  4a.  The  flow  for  this  roughness 
topology  was  expected  to  follow  the  canonical  mean  and  turbulent  patterns  observed  by  Goddard 
(1959)  and  Latin  and  Bowersox  (2000),  respectively. 
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0.8  mm 


(a)  Square  Pattern 


(b)  Diamond  Pattern 


Fig.  3  Roughness  Topology. 


(a)  Square  pattern  (b)  Diamond  Pattern 

Fig.  4  Sketch  of  the  square  (side  view)  and  diamond  (top  view)  roughness  element  flow  patterns. 


The  diamond  roughness  topology  was  seleeted  so  that  the  height  and  frontal  bloekage  elosely 
matehed  the  square  roughness  elements.  The  diamond  roughness,  however,  indueed  signifieant 
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localized  shock  and  expansion  flow  distortions  by  way  of  attached  leading  edge  oblique  shoek 
waves,  followed  by  Prandtl-Meyer  expansions.  The  diamond  elements  were  9.0  mm  long,  1.6 
mm  wide  and  0.8  mm  tall.  The  major  axis  of  the  roughness  element  was  aligned  with  the  flow 
direetion,  and  the  leading  edge  had  a  10°  half  angle.  A  sehematic  of  the  element  shape  is  given  in 
Fig.  3b,  while  the  expected  flow  structure  is  illustrated  in  Fig.  4b.  As  indicated  in  Fig.  4b,  the 
flow  was  expeeted  to  wrap  around  the  elements,  thus  exposing  more  of  the  topology  to  the 
supersonie  cross-flow.  If  the  entire  element  was  exposed  to  the  flow,  then  the  average  loeal 
Maeh  number  at  the  element  height  would  have  been  1.3.  Sinee  the  boundary  layer  was 
adiabatic,  the  shoek  turbulenee  interaction  was  expected  to  be  similar  to  ease  IB  in  Mahesh  et  al. 
(1997). 


2.3  Pressure  gradients 

Two  pressure  gradient  models  were  designed  to  provide  signifieant  pressure 
gradient/streamline  curvature  effects.  Expansion  and  compression  wave  theories  were  used  to 
estimate  the  values  of  the  distortions  and  impulse  parameters  for  the  global  pressure  gradients. 
The  Weak  Favorable  Pressure  Gradient  (WFPG)  model  was  designed  to  provide  a  small 
perturbation,  and  the  Strong  Favorable  Pressure  Gradient  (SFPG)  model  produced  a  more 
pronouneed  effect,  where  the  imposed  global  gradient  was  of  similar  magnitude  but  opposite 
sign  to  the  adverse  gradient  on  the  front  half  of  the  diamond  surfaee  element.  The  eombined 
loeal  and  global  gradients  were  ehosen  to  help  isolate  the  shoek  turbulence  interaction  effects. 
The  pressure  gradient  models  were  similar  to  those  used  by  Luker  et  al.  (2000).  The  eurved  walls 
were  maehined  to  match  the  following  polynomial:  y  =  A(x-xi)^  +  5(x-xi)^,  where  xi  is  the  axial 
loeation  where  the  eurvature  begins.  The  beginning  and  end  loeations  for  the  wall  curvature. 
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along  with  the  polynomial  coefficients,  are  listed  in  Table  1.  A  schematic,  with  the  coordinate 
system,  is  given  in  Fig.  5.  A  flat  wall  was  used  to  approximate  a  Zero  Pressure  Gradient  (ZPG) 
flow.  In  addition  to  the  favorable  pressure  gradient  region,  measurements  were  also  acquired  just 
downstream  of  the  compression  region  for  the  strong  favorable  pressure  gradient  model,  which 
will  be  referred  to  as  the  Combined  Pressure  Gradient  (CPG)  region. 


Table  1:  Model  geometry  and  pressure  gradient  parameters 


X;  (cm) 

X2  (cm) 

A 

B 

WFPG 

24.54 

63.54 

6.423  X  10'^ 

-8.757  X  10'^ 

SFPG 

24.54 

39.64 

1.107  X  10"^ 

-2.507  X  10'^ 

2.4  Coordinate  system  and  measurement  locations 

A  wind  axes  Cartesian  coordinate  system  was  adopted,  where  the  origin  was  located  on  the 
floor  at  the  center  of  the  test  section  entrance.  The  coordinates  in  Fig.  5  correspond  to  the 
adopted  system.  The  roughness  element  plates  started  at  x  =  0.  Measurements  were  centered  on 
three  locations;  x  =  15.9  cm,  29.8  cm  and  41.3  cm  (as  shown  in  Fig.  5).  For  orientation  purposes, 
the  leading  edge  of  the  18*  row  of  diamond  elements  for  all  of  the  models  was  located  at  x  = 
15.5  cm.  This  location  corresponded  to  the  49*  row  of  square  roughness  elements. 


2.5  Schlieren  photography 

A  Z-type  2-mirror  schlieren  system  was  used.  The  knife-edge  was  aligned  for  density  gradients 
normal  to  the  wall.  The  schlieren  image  was  recorded  by  a  Nikon  D50  SLR  digital  camera. 


X  (cm) 

Fig.  5  Sketch  of  streamline  curvature  wall  geometry. 
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2.6  Pressure  sensitive  paint 

Surface  pressure  distributions  were  doeumented  using  Pressure  Sensitive  Paint  (PSP)  [Liu 
and  Sullivan  (2005)].  The  paint  used  in  the  present  tests  was  a  UNC405  produced  by  Innovative 
Seientifie  Solutions  Ine.  (ISSI).  The  binary  paint  was  designed  for  a  pressure  range  of  10  and  150 
kPa,  whieh  eovered  the  surfaee  pressures  in  this  study.  Six  coats,  each  approximately  10  pm 
thiek,  were  sprayed  onto  each  test  seetion  model.  The  models  were  eured  with  a  heat  gun  for 
approximately  15  minutes  after  the  final  coat  was  applied.  An  ISSI  Ine.  Lm  series  LED  array 
light  souree  was  used  to  exeite  the  PSP  and  a  PCO  1600  camera  was  used  to  image  the 
luminescenee  over  5.08  cm  by  3.81  cm  area.  The  1600  by  1200  eamera  resolution  was  redueed 
to  800  by  600  via  pixel  binning  for  enhaneed  signal  to  noise  ratio. 

OMS  Lite  version  1.1  by  ISSI  was  the  software  paekage  used  to  proeess  the  PSP  images.  The 
software  was  eapable  of  resolving  errors  due  to  photo-deteetor  noise,  non-uniform  image 
illumination  and  model  deformation  by  using  a  two  dimensional  displaeement  evaluation  which 
included  image  registering.  For  eaeh  run,  four  images  were  acquired:  a  baekground  image 
without  ambient  illumination,  an  initial  wind-off  image  taken  just  before  the  tunnel  was  started 
and  with  the  LED  array  on,  a  wind-on  image  with  the  tunnel  running  long  enough  to  reaeh 
thermal  equilibrium,  and  a  final  wind-off  image  taken  immediately  after  the  tunnel  was  stopped. 
To  account  for  model  distortions,  8  marker  dots  were  applied  to  the  surfaee  with  a  felt  tip  pen  at 
known  loeations.  The  loeation  of  each  marker  dot  was  entered  into  the  OMS  software  and  image 
registration  was  aceounted  for  using  an  8-step  Quantum  Pixel  Energy  Distribution  (QPED) 
algorithm  as  deseribed  by  Navarra  et  al.  (2001)  with  3  iterations  and  an  initial  box  size  of  16  by 
16  pixels.  Due  to  the  paint’s  thermal  and  pressure  dependenee,  the  PSP  was  ealibrated  by  a 
seeond-order  least-squares  approximation  with  known  temperatures  and  pressures.  Finally  a  4- 
pixel  box  smoothing  filter  was  applied  to  sharpen  the  linear  interpolation  of  PSP  values. 
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Comparisons  were  made  on  the  smooth  model,  and  the  PSP  values  were  within  1.0%  of  the 
pressure  tap  data.  However,  due  to  the  small  size  of  the  roughness  elements,  it  was  not  possible 
to  resolve  the  pressure  on  the  surfaee  of  the  elements  or  within  the  ehannels  using  statie  pressure 
ports.  The  measurement  uneertainty  analyses  for  the  PSP,  as  well  as  the  remaining  diagnostic 
methods,  are  summarized  in  Section  2.9. 

2.7  Pitot  pressure 

The  pitot  pressure  was  measured  with  an  Omega  brand  PX303  pressure  transducer  (0  -  690 
kPa,  and  +0.25%  full-scale  accuracy).  The  pressure  probe  had  a  1.6  mm  diameter  tube  with  an 
inner  diameter  of  0.8  mm.  The  end  of  the  tube  was  crimped  so  that  the  entrance  height  was  0.2 
mm.  The  probe  was  traversed  in  the  floor-normal  direction  by  a  Velmex  model  MA2508P40 
linear  translation  stage  that  contained  a  turn  screw  that  had  15.75  turns  per  centimeter.  The 
translation  stage  was  moved  by  a  Vexta  PK2608-02A  400-step  stepper  motor.  The  stepper  motor 
was  operated  by  an  SI 00  3 -channel  stepper  motor  controller  which  was  controlled  by  a  PC 
running  the  SI 00  SMC  software.  The  probe  was  translated  at  1.9  cm  per  minute.  The  position  of 
the  probe  was  monitored  with  a  P3  America  Inc.  ISDG  linear  variable  displacement  transducer 
(0  -  20  mm  and  +0.5%  full-scale  accuracy).  Probe  pressures  were  sampled  at  a  rate  of  1.0  kHz, 
and  each  data  point  was  the  result  of  a  100-point  average. 

2.8  Particle  image  velocimetry 

Particle  image  velocimetry  was  used  to  measure  velocity  fields  across  the  boundary  layer 
with  sufficient  spatial  resolution  to  resolve  the  flow  over  the  roughness  elements.  For  the  present 
tests,  the  tunnel  was  seeded  by  a  TSI  Six-jet  Atomizer  Model  9306  with  Bertolli  extra  light  olive 
oil.  The  documented  mean  particle  diameter  was  0.6  pm,  and  the  3  dB  frequency  response  of  the 
particles  was  documented  at  60  -  200  kHz  [Bowersox  (1996)].  Tracer  particles  were  illuminated 
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via  Mie  scattering  from  laser  radiation  supplied  by  a  dual  port/dual  head  New-Wave  Solo  120XT 
PIV  laser.  Each  laser  head  had  a  maximum  energy  output  of  120  mJ  at  532  nm,  4  ns  pulse  width 
and  a  maximum  rep-rate  of  15  Hz.  To  minimize  laser  energy  losses  during  beam  collimation,  the 
beam  out  of  the  second  port  was  cross-polarized  using  a  crystal  quartz  half-wave  plate,  and 
collimated  using  a  high-energy  polarizing  cube  beam-splitter  that  provided  efficient  narrowband 
polarization.  A  BK7  spherical  lens  {f=  500  mm)  was  used  to  create  a  nominally  500  |am  thick 
laser  sheet,  above  the  test  section  floor.  The  sheet  width  was  around  3.0  cm,  and  was  formed 
using  a  BK7  Plano-concave  cylindrical  lens  (/’=  -  75  mm).  The  laser  sheet  was  located  along  the 
lateral  centerline  of  the  wind  tunnel.  The  camera  used  was  a  Cooke  Corp.  PCO  1600  interline- 
transfer  CCD  camera  equipped  with  a  Nikon  60  mm  lens.  The  camera  trigger,  laser  Q-switch  and 
laser  flashlamp  were  all  controlled  by  a  Quantum  Composers  Model  9618  pulse  generator.  The 
camera  frame  grabbing  software  was  Camware  version  2.13. 

Velocity  fields  were  created  using  the  ISSl  brand  dPlV  32-bit  Analysis  Code.  For  the  present 
measurements,  a  three-step  adaptive  cross  correlation  was  used,  with  successive  rectangular 
interrogation  sizes  of  128x64,  64x32,  and  32x16  pixels  and  a  50%  overlap.  In  order  to  enhance 
the  intensity  of  correlation  peaks  relative  to  random  noise,  a  four-map  correlation  multiplication 
process  filter  was  used.  A  consistency  post-processing  filter  was  used  to  improve  the  adaptive 
correlation  calculation.  First  and  second  order  turbulent  statistics  were  created  using  an  in-house 
code  that  ensemble  averaged  the  velocity  vector  fields.  For  the  present  study,  5000  image  pairs 
were  acquired  at  each  measurement  location.  The  measurement  grid  consisted  of  71  by  34  data 
points,  where  the  axial  and  transverse  separations  were  0.71  mm  and  0.35  mm,  respectively.  This 
resolution  was  sufficient  to  compute  gradients  of  the  mean  and  turbulent  flow  properties  across 
the  roughness  elements. 
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2.9  Measurement  Uncertainties 


Measurement  uneertainties  for  the  present  study  are  summarized  in  Table  2,  and  were 
accumulated  with  a  Euclidean  {Li)  norm.  Uncertainties  in  stagnation  conditions  included  probe 
location  uncertainty  and  run-to-run  repeatability.  The  run-to-run  repeatability  of  the  PSP  was  3.3 
kPa.  However,  the  uncertainty  between  different  measurement  points  for  the  same  run  was 
considerably  smaller  (~200  Pa).  The  axial  position  uncertainty  was  taken  as  half  the  smallest 
division  on  the  calibration  scale,  and  the  pitot  probe  transverse  probe  position  uncertainty  was 
taken  as  approximately  half  the  probe  tip  height.  Boundary  layer  thickness  uncertainty  accounted 
for  the  uncertainties  in  position  and  velocity,  where  the  sensitivity  to  the  velocity  uncertainty  was 
estimated  with  a  1/7*  power  law.  The  uncertainty  analysis  of  the  PIV  data  accounted  for 
correlation  mapping  error  and  the  conversion  error  from  the  physical  length  scale  to  the 
appropriate  number  of  camera  pixels.  The  axial  position  uncertainty  was  the  same  as  that  of  the 
pitot  probe.  The  transverse  position  uncertainty  was  taken  as  that  of  the  wall  position,  which  was 
known  to  within  2  pixels  (~0.1  mm)  on  the  camera  array.  To  estimate  the  uncertainty  in  the 
statistical  quantities,  a  95%  confidence  interval  was  used  [Benedict  and  Gould  (1996)],  with  the 
variance  determined  assuming  a  normal  distribution  of  5000  instantaneous  velocity  vector  fields. 
The  normalized  uncertainty  analysis  was  only  applicable  to  the  axial  velocity  because  magnitude 
of  the  transverse  velocity  was  comparatively  very  small.  Hence,  the  magnitude  of  the  transverse 
velocity  uncertainty  was  taken  to  be  the  same  as  that  of  the  axial  velocity.  The  length  scale  and 
axial  velocity  uncertainties  for  the  PIV  and  pitot  pressure  probe  were  similar.  The  tabulated 
fluctuating  velocity  correlation  uncertainty  is  for  the  boundary  layer  outer-region  only;  near  the 
wall  the  uncertainties  were  considerably  higher,  e.g.,  -35%  ditylS-  0.1.  The  magnitude  of  the 
velocity  gradient  uncertainty  was  a  combination  of  the  uncertainties  in  spatial  position  and  mean 
velocity  measurements.  The  spatial  difference  uncertainty  only  accounted  for  the  conversion 
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error  from  the  physical  length  scale  to  the  PIV  image  pixels  as  the  offset  uncertainties  applied  to 
the  entire  image.  The  mean  velocity  uncertainty  included  the  value  listed  in  Table  2.  Added  to 
the  spatial  difference  uncertainty  for  y-component  derivative  was  an  estimate  of  the  error  due  to 
the  curvature  of  the  velocity  profile  within  PIV  interrogation  window.  The  resulting  uncertainty 
estimates  are  listed  in  Table  2.  The  resulting  uncertainty  for  the  production  was  estimated  at 
approximately  20%. 


Table  2:  Uncertainty  Estimates 


Variable 

Error 

poi  ?Toi 

3.0  [kPa],  16  [K] 

PSP 

3.3  [kPa]" 

X,  y 

0.4,  0.1  [mm]*’ 

s 

0.7  [mm]’’ 

u 

3.8  [m/s]’’ 

Ueff 

6.0%’’ 

u'u' ,  v'v' 

4.0% 

m'v' 

8.0% 

du  !  dx  ,du  !  dy 

3.0%,  8.0% 

P^IP^  Pyylp^P^lp 

20.0% 

^Refers  to  mn-to-run  repeatability.  Spatial  uneertainty  from  the  same  run  was  mueh  smaller  (~200  Pa) 

^Applies  to  both  pitot  probe  and  PIV 


2.10  Flow  Field  Three-Dimensionality 

Elimination  of  laser  reflections  off  of  the  roughness  elements  was  a  key  challenge  associated 
with  the  present  study,  where  data  were  acquired  along  the  centerline  of  the  tunnel.  For  the 
smooth  and  square  roughness  pattern,  it  was  reasonable  to  assume  that  the  flow  statistics  were 
nominally  two-dimensional,  and  thus  the  laser  sheet  was  located  such  that  wall  reflections  were 
minimized.  This  assumption  was  questionable  for  the  diamond  roughness  elements,  as  localized 
secondary  flows  were  expected.  To  estimate  the  effects  of  three-dimensionality  on  the  diamond 
element  flow,  additional  PIV  measurements  were  taken  across  the  span.  The  first  measurements 
were  taken  along  the  center  of  the  lateral  span  of  a  roughness  element,  while  a  second  data  set 
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was  acquired  at  a  position  shifted  0.4  mm  in  the  z-direction,  which  was  near  the  midpoint 
between  the  element  center  and  element  edge.  The  shifted  data  were  subtracted  from  the 
centerline  contours.  The  transverse  velocity  difference  contour  is  displayed  in  Fig.  6.  The  dark 
vertical  structure  located  near  15.75  and  the  bright  axial  structure  neary/t^e/=  0.15  were 

the  results  of  laser  reflections.  Away  from  these  reflections,  the  maximum  difference  was 
nominally  3.0  m/s,  which  was  less  than  the  3.8  m/s  measurement  uncertainty  listed  in  Table  2. 
The  axial  velocity  had  a  similar  difference.  Therefore,  for  the  present  work,  the  relatively  small 
three-dimensional  effects  were  neglected  and  the  laser  sheets  were  positioned  to  minimize  laser 
reflections. 


0  0 . .  .  .  1  . 

15.0  15.5  16.0  16.5 


0.025 

0.020 

0.015 

0.010 

0.005 

O.QQO 

-0.QQ5 

-0.010 

-0.015 

-0.020 

-0.025 


X 


IS, 


ref 


Fig.  6  Diamond  roughness  element  velocity  difference:  (v^.^  ~^side)l^ref  * 


3.  Overall  Flow  Structure 


3.1  Free  stream  and  Surface  Properties 

The  measured  flow  and  surface  properties  for  the  present  study  are  summarized  in  Table  3. 
The  freestream  Mach  numbers  and  boundary  layer  thicknesses  were  based  on  pitot  probe 
surveys.  The  wall  density  and  temperature  were  computed  assuming  an  adiabatic  wall.  The 
friction  velocity  was  estimated  by  assuming  that  the  Townsend  (1976)  trend,  with  Morkovin 
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(1961)  scaling,  shown  in  Fig.  lb,  held  for  the  zero  pressure  gradient  smooth  and  J-type 
roughness  flows.  The  present  Reynolds  shear  stress  data  results  for  the  smooth,  square  and 
diamond  model  eases  are  ineluded  in  Fig.  lb.  In  this  plot,  the  wall  shear  was  seleeted  to  best  fit 
the  trend  near  the  wall.  It  was  assumed  that  the  near  wall  behavior,  where  t  ,  +  ,  =  T„, ,  was 

most  relevant.  The  Reynolds  shear  stress  data  were  averaged  axially  for  a  distanee  of  1.64-e/to 
minimize  the  effects  of  the  loeal  roughness  generated  pressure  gradients  on  the  wall  shear  stress 
estimates.  As  indicated  in  Fig.  lb,  the  overall  shapes  of  the  present  smooth  and  square  roughness 
profiles  followed  the  Klebanoff  (1959)  trend,  as  did  the  smooth  wall  supersonic  data  of  Luker  et 
al  (2000).  However,  the  diamond  stresses  were  30%  below  the  trend  in  the  outer-region  of  the 
boundary  layer.  This  difference,  like  those  described  in  Latin  and  Bowersox  (2000),  was  due  to 
the  influence  of  the  loeal  shocks  and  expansions  on  the  turbulent  shear  stresses;  these  phenomena 
are  described  in  more  detail  below.  The  resulting  frietion  velocities  are  listed  in  Table  3a. 


Table  3a:  Flow  Conditions  for  Zero  Pressure  Gradient  Models 


Model 

^oo 

s 

mm 

Ur 

m/s 

Pw 

kg/m^ 

T 

j-  w 

K 

kjk 

K 

X 

=  15.9  cm 

ZPG  Smooth  Plate 

2.85 

11 

23 

0.28 

294 

0 

0 

ZPG  Square  Roughness 

2.82 

12 

35 

0.29 

294 

0.73 

440 

ZPG  Diamond  Roughness 

2.84 

11 

26 

0.28 

294 

0.07 

23 

21 


Table  3b:  Flow  Conditions  for  Non-Zero  Pressure  Gradient  Models 


Model 

d 

mm 

Pw 

kg/m^ 

T 

j-  w 

K 

X  =  29.8  cm 

WFPG  Smooth  Plate 

2.87 

13 

0.27 

294 

WFPG  Square  Roughness 

2.84 

14 

0.28 

294 

WFPG  Diamond  Roughness 

2.87 

14 

0.27 

294 

SFPG  Smooth  Plate 

3.16 

16 

0.18 

292 

SFPG  Square  Roughness 

3.12 

17 

0.19 

293 

SFPG  Diamond  Roughness 

3.14 

16 

0.18 

292 

X  =  41.3  cm 

WFPG  Smooth  Plate 

2.97 

16 

0.23 

293 

WFPG  Square  Roughness 

2.90 

17 

0.26 

294 

WFPG  Diamond  Roughness 

2.94 

17 

0.24 

293 

CPG  Smooth  Plate 

2.77 

19 

0.31 

294 

CPG  Square  Roughness 

2.81 

23 

0.30 

294 

CPG  Diamond  Roughness 

2.81 

20 

0.30 

294 

3.2  Schlieren  photography 

The  estimated  wall  shear  stress  for  the  square  roughness  elements  was  nearly  twiee  that  for 
the  diamond  elements.  Following  Goddard’s  (1959)  arguments,  a  lower  wall  shear  stress  was 
expeeted  for  the  more  streamlined  diamond  elements.  The  eorresponding  equivalent  sand  grain 
heights  were  computed,  as  described  in  Section  1.1,  and  included  in  Table  3a.  Wall  shear  stresses 
were  not  estimated  for  the  cases  with  imposed  global  distortions  (Table  3b),  as  the  presence  of 
the  pressure  gradient  would  have  introduced  significant  uncertainty. 

Schlieren  photographs  were  acquired  to  provide  an  overall  view  of  the  flow  structure. 
Because  of  optical  access  limitations,  images  were  restricted  to  the  second  station  (x  =  29.8  cm). 
Examples  for  the  strong  favorable  pressure  gradient  model  are  given  in  Figs.  7a  and  7b.  The 
local  roughness  induced  flow  structure  for  the  weak  and  zero  pressure  gradient  models  followed 
the  patterns  described  below.  The  flow  direction  was  from  left  to  right  and  the  surface  roughness 
was  located  along  the  lower  surface.  Surface  boundary  layers  in  Fig.  7a  are  visible  as  the  dark 
and  light  bands  along  the  floor  and  ceiling,  respectively.  The  dark  line  angled  down  to  the  right 
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emanating  from  about  two-thirds  the  way  up  from  the  bottom  left  comer  was  a  weak  wave  which 
formed  at  the  intersection  of  nozzle  exit  and  test  section  floor,  and  was  too  weak  to  be  detected 
by  the  quantitative  measurements.  In  addition,  faint  waves  were  observed  originating  from  the 
square  roughness  elements.  The  diamond  element  roughness  model  (Fig.  7b)  produced  a 
significantly  different  wave  stmcture  than  was  seen  for  the  square  roughness  model.  The  shock 
waves  originating  from  the  roughness  elements  are  clearly  discernible  as  the  bright  curved  lines 
that  moved  up  and  to  the  right.  These  waves  were  significantly  stronger  than  those  for  the  square 
roughness  model.  The  flow  expansion  for  this  model  was  also  visible  as  shock  waves  spread 
apart  as  the  flow  moved  around  the  wall  curvature.  The  PIV  contours  (Fig.  10b  and  d)  for  the 
diamond  element  models  displayed  a  wave  pattern,  which  corresponded  to  the  wave  stmcture 
shown  in  Fig.  7b;  these  waves  were  also  present  in  the  pitot  surveys  (not  shown).  These  wave 
patterns  were  not  present  in  the  corresponding  square  roughness  data.  These  images  confirm  the 
presence  of  roughness  element  generated  gas  dynamic  waves. 


(a)  Square  Roughness  Topology  (b)  Diamond  Roughness  Topology 

Fig.  7  Schlieren  photograph  at  2"**  station  for  strong  favorable  pressure  gradient  roughness  models 

3.3  Pressure  sensitive  paint 
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The  PSP  data  were  used  primarily  to  examine  the  loeal  pressure  fields  generated  by  the 
roughness  elements.  Example  surface  pressure  distributions  for  the  square  and  diamond 
roughness  patterns  are  presented  in  Figs.  8a  for  the  zero  and  strong  favorable  pressure  gradient 
models,  with  corresponding  line  traces  for  the  zero  pressure  gradient  models  plotted  in  Fig.  8b. 
For  the  square  roughness  topology,  a  modest  5%  increase  in  surface  pressure  was  observed 
within  the  channels  between  the  roughness  element  surfaces,  while  the  diamond  roughness 
model  PSP  image  revealed  strong  alternating  adverse  and  favorable  pressure  gradients  which 
were  not  restricted  to  the  interface  between  the  roughness  edge  and  the  channels.  These  data 
confirmed  the  expected  flow  topology  shown  in  Fig.  4b,  where  a  sketch  is  overlaid  on  the  PSP 
image.  On  the  front  half  of  the  roughness  element,  there  was  an  increase  in  the  pressure  due  to 
the  oblique  shock  waves.  The  pressure  then  decreased  on  the  back  half  as  a  result  of  the 
expansions.  The  high-  and  low-pressure  regions  from  the  individual  elements  overlapped 
adjacent  elements  producing  pressure  gradients  that  exhibited  a  strong  periodic  axial  variation 
(±10  to  18%  of  the  mean  wall  pressure).  The  amplitude  of  the  pressure  rise  was  found  to  depend 
on  the  position  relative  to  the  roughness  elements.  The  two  larger  peaks  at  x/5ref=  29.5  and  30.7 
in  Fig.  8b  corresponded  to  the  regions  between  the  leading  edges  of  two  adjacent  elements,  while 
the  smaller  one  located  at  x/5ref  =  30.0  corresponded  to  the  pressure  directly  on  a  roughness 
element.  The  strong  favorable  pressure  gradient  models  had  similar  localized  flow  structure 
around  the  roughness  elements,  but  also  exhibited  the  expected  steep  decrease  in  surface  pressure 
in  the  axial  direction.  Comparing  the  left  and  right  images  in  Fig.  8a,  the  noticeable  change  in 
grayscale  indicates  the  global  pressure  reduction  was  substantial. 
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(a)  ZPG  (left)  and  SFPG  (right)  Roughness  Pressure  Contours  (b)  Line  Profiles  (ZPG  case) 

Fig.  8  Pressure  Sensitive  Paint  results 

3.4  Zero  Pressure  Gradient  Inner  and  Outer  Law  Scaled  Velocity  Plots 

Law-of-the-wall  and  defect  law  plots  were  constructed  for  the  zero  pressure  gradient  models. 
The  velocity  profiles  were  extracted  from  the  PIV  data,  and  computed  from  the  pitot  data.  The 
Mach  number  was  estimated  from  the  pitot  data  using  the  Rayleigh-Pitot  formula,  where  the 
static  pressure,  which  was  assumed  constant,  was  computed  based  on  the  measured  freestream 
total  pressure  ratio.  Below  the  sonic  point,  the  Mach  number  was  computed  from  the  ratio  of  the 
local  total  and  static  pressure  using  the  isentropic  relations.  The  velocity  was  then  estimated 
assuming  a  constant  total  temperature  across  the  boundary  layer.  Latin  and  Bowersox  (2000) 
demonstrated  the  feasibility  of  this  approach  for  smooth  and  rough  supersonic  boundary  layers 
under  similar  flow  conditions.  For  the  present  measurements,  the  planar  PIV  data  were  averaged 
in  the  axial  direction  over  a  range  of  ~1.6  <^e/,  corresponding  to  10  square  roughness  elements,  to 
minimize  the  effects  of  the  local  roughness  topology.  Smooth  and  square  roughness  data  were 
minimally  affected  by  the  averaging  process,  where  the  scatter  was  within  the  measurement 
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uncertainty.  However,  for  the  diamond  case,  the  velocity  and  near  wall  Reynolds  shear  stress 
varied  significantly  across  the  roughness  element.  Thus,  the  diamond  model  velocity  data  are 
included  in  this  section  for  qualitative  comparisons  only.  In  addition,  the  data  was  not  spanwise 
averaged,  further  lowering  the  validity  of  the  diamond  model  averages. 

Plots  with  inner  and  outer  scaling  are  displayed  in  Figs.  9a  and  9b,  respectively.  The  van 
Driest  II  compressibility  scaling  was  also  employed,  where  1  = 

I uj)  — B^]! {B^  I B^)  and  the  5,  are  included  in  the  nomenclature 

section.  Skin  friction  was  estimated  from  the  PIV  turbulence  shear  stress  data  as  described  in 
Section  3.1.  Smooth  model  PIV  and  pitot  data  exhibited  good  agreement  with  each  other  and  the 
law-of-the-wall  correlation  (Fig.  9a)  given  by  =  K~^\n,y^  +B  +  211  K~^  sin^  {tt  !  1- y !  S) , 

where  the  constants  K,  B  and  77  were  set  to  0.4,  5.5  and  1.0,  respectively.  The  square  roughness 
model  produced  the  expected  downward  shift  of  the  law-of-the-wall  for  both  the  PIV  and  pitot 
data,  again  showing  very  good  agreement  with  each  other.  The  downward  shift  was  «  12.2. 

The  diamond  model  PIV  data  are  also  plotted  in  Fig.  9a,  and  the  downward  shift  was  estimated 
at  4.8.  This  estimate  was  considered  qualitative  for  the  reasons  listed  in  the  previous  paragraph 
and  because  only  three  data  points  were  present  in  the  logarithmic  layer.  The  intercepts  C  of 
hr y  / A: -I- C  were  9.3  and  15.0  for  the  square  and  diamond  roughness  models, 

respectively.  The  equivalent  sand  grain  heights  were  then  computed  following  Schlichting 
(1979),  and  the  results  are  summarized  in  Table  3a.  These  data  are  also  plotted  in  Fig.  la;  labeled 
Present  (M  =  2.8).  The  smooth,  axially  averaged  diamond  and  square  roughness  velocity  data  all 
collapsed  onto  the  defect-law  plot  (w^^  ~^eff  2nx'“'[l-sin^(.;r/  l-ylS)-\, 

when  van  Driest  II  scaling  was  used  (Fig.  9b). 
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(a)  Inner  Scaling  (b)  Outer  Scaling 

Fig.  9  Velocity  profiles  for  the  zero  pressure  gradient  smooth,  square  and  diamond  roughness  models. 


3.5  Rough  Wall  Velocity  Contours 

Mean  velocity  contours  for  the  zero  pressure  gradient  square  and  diamond  roughness  model 
are  shown  in  Fig.  10.  The  data  in  Fig.  10  were  representative  of  all  of  the  zero  pressure  gradient 
locations.  The  velocity  components  and  length  scales  were  normalized  by  the  reference  velocity 
(625  m/s)  and  reference  boundary  layer  thickness  (1.0  cm),  respectively. 

The  zero  pressure  gradient  square  roughness  u  velocity  contours  (Fig.  10a)  were  uniform  in 
the  axial  direction  and  decreased  in  the  transverse  direction  in  the  expected  manner,  while 
corresponding  v  -components  were  very  small.  On  the  other  hand,  the  u  -component  contour 
plots  for  the  diamond  roughness  exhibited  a  wavy  pattern  throughout  the  boundary  layer  and  into 
the  ffeestream.  The  pattern  corresponded  to  shock  and  expansion  waves  generated  by  the 
roughness  elements  as  shown  in  Fig.  7b.  The  low  velocity  regions  near  the  wall  at  x/<^e/=15.7 
and  16.6  were  near  the  leading  edge  of  a  row  of  diamond  elements.  In  addition,  for  the  v  - 
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component,  a  series  of  alternating  negative  and  positive  bands  were  observed,  where  the  values 
varied  between  0.02  and  -0.01 5  [4-e/. 


(c)  Square  Roughness  (  V  /  U^^j- )  (d)  Diamond  Roughness  ( V  /  ) 

Fig.  10  Zero  Pressure  Gradient  Rough  Wail  Velocity  (PIV)  Contours  (Location  1) 


The  strong  favorable  pressure  gradient  model  transverse  veloeity  results  for  the  smooth, 
square  and  diamond  roughness  element  plates  are  given  in  Fig.  11.  The  weak  favorable  pressure 
gradient  results  were  very  similar  to  the  zero  pressure  gradient  results,  and  thus  the  contours  are 
not  reprodueed  here.  However,  the  trends  are  ineluded  in  the  subsequent  diseussion.  The  mean 
transverse  veloeity  eomponent  ( v  )  at  the  second  station  for  the  weak  favorable  pressure  gradient 
model  was  nominally  -O.OlUref  for  both  the  smooth  and  square  surface  roughness  models,  with 
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similar  trends  exhibited  at  the  third  station.  The  only  difference  of  note  between  the  second  and 
third  stations  was  the  larger  turning  angle,  where  v  =  -0.05  f/re/  for  the  smooth  and  -0.04 17^/  for 
the  square  roughness  cases,  which  are  in  line  with  the  steeper  wall  angle  (see  Fig.  5).  The  strong 
favorable  pressure  gradient  smooth  wall  transverse  velocities  varied  between  -0.05  and  -0.15L7e/ 
over  the  PIV  interrogation  window  at  Station  2.  The  corresponding  square  roughness  model 
boundary  layer  was  thicker  and  the  transverse  velocity  range  was  slightly  smaller  at  -0.05  to  - 
QAlUref-  In  contrast,  by  the  combined  pressure  gradient  region  (Fig.  11b)  the  flow  had  resumed 
in  the  axial  direction,  as  indicated  by  the  nominally  zero  transverse  velocity  component.  It  should 
be  noted  however,  that  in  the  upper  left  comer,  which  was  located  upstream  of  the  recompression 
shock,  the  transverse  velocity  was  still  moving  towards  the  floor.  The  axial  velocity  contours 
(not  shown)  were  qualitatively  similar  for  all  three  models  and  similar  in  stracture  to  the 
upstream  zero  pressure  gradient  contours  shown  in  Figs.  10a  and  b.  For  the  square  model,  a 
region  of  relatively  low  axial  momentum  was  present  in  the  lower  40%  of  the  boundary  layer. 
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Fig.  11  SFPG  Model  Transverse  Velocity  (PIV)  Contours 
For  the  diamond  roughness  cases  with  favorable  pressure  gradients,  the  u  -component 

contour  plots  exhibited  a  wavy  pattern  similar  to  that  observed  in  Fig.  10  for  the  zero  pressure 

gradient  models.  Also,  low-momentum  regions  were  visible  near  the  leading  edges  of  two  rows 

elements  at  x/(^e/=29.2  and  30.2  of  the  weak  favorable  pressure  gradient  and  strong  favorable 

pressure  gradient  models.  The  v  contour  for  the  second  weak  favorable  pressure  gradient  region 

(not  shown)  displayed  the  same  periodic  pattern  seen  in  Fig.  11;  however,  in  the  expansion 

region,  the  velocity  range  was  from  0  to  -O.Odf^e/.  If  v  due  to  streamline  curvature,  estimated  at 

-0.02Ur^,  based  on  the  smooth  and  square  element  roughness  results,  is  subtracted  from  the 

weak  favorable  pressure  gradient  diamond  roughness  values  element  values  of  v  ,  a  wave  pattern 

resulted  with  a  relative  variation  due  to  the  roughness  elements  that  was  similar  to  the  zero 

pressure  gradient  case.  Again,  aside  from  the  increase  in  flow  turning  angle,  the  trends  at  the 
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third  region  were  qualitatively  similar.  For  the  strong  favorable  pressure  gradient  model  at 
Loeation  2,  the  variation  observed  between  v  ranged  between  -0.07  and  -Q.\5Uref-  The  greater 
peak-to-peak  variation,  relative  to  the  zero  pressure  gradient  case,  indicated  enhanced  locally 
generated  strain-rates  with  the  globally  imposed  pressure  gradient.  The  enhanced  strain-rates  will 
be  discussed  in  more  detail  in  the  following  section.  For  the  combined  pressure  gradient  region, 
the  periodic  nature  of  the  transverse  velocity  component  was  noticeably  absent  for  the  diamond 
model. 

3.6  Section  Summary 

Diamond,  A:-type,  roughness  altered  the  flow  in  a  manner  that  was  significantly  different  than 
the  square,  t/-type,  roughness  of  comparable  height  and  frontal  blockage.  Specifically,  the 
diamond  element  shock  and  expansion  waves  had  a  measureable  effect  on  the  local  surface 
pressure  and  mean  flow  field.  The  flow  structure  for  the  zero  pressure  gradient  square  roughness 
topology  followed  the  sand  grain  behavior  described  by  Goddard  (1959).  When  spatially 
averaged,  the  diamond  element  data  also  tended  to  collapse  onto  the  law  of  the  wall  and  defect 
law  plots,  with  van  Driest  II  (1951)  compressibility  scaling.  The  diamond  element  model 
collapse  is  not  considered  definitive  as  the  averages  were  insufficient.  Nonetheless,  this  tendency 
suggests  that  the  average  behavior  of  the  diamond  element  flow,  even  with  the  strong  wave 
structure,  is  amenable  to  the  engineering  approximations  described  in  Section  1.2.  The  favorable 
pressure  gradients  strengthened  the  secondary  flow  structure  for  the  diamond  roughness 
elements.  However,  in  the  combined  pressure  gradient  region,  the  sonic  point  evidently  moved 
above  the  roughness  height,  and  wave  structure  was  absent.  Thus,  in  this  region,  the  A:-type 
diamond  elements  were  expected  to  behave  more  like  their  incompressible  counterparts. 
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4.  Response  of  the  Reynolds  Stresses  to  the  Fluid  Element  Deformations 


Turbulence  stresses,  in  the  measurement  coordinate  system,  were  computed  from  the  PIV 
data  for  all  of  the  cases  studied.  In  this  section,  the  response  of  the  Reynolds  shear  stresses  to  the 
mean  fluid  element  deformations  is  characterized.  The  normal  stresses,  in  principal  axes,  are 
discussed  in  Section  5.  For  the  present  discussion,  the  canonical  flows  were  defined  as  the  cases 
where  the  boundary  layer  flow  structure  depends  on  the  surface  condition  and/or  pressure 
gradient  in  a  global  manner.  Conversely,  the  non-canonical  flows  were  defined  as  the  cases 
where  the  local  boundary  flow  structure  was  significantly  altered  by  the  wall  topology.  The 
canonical  flows  include  all  of  the  smooth  wall  cases,  all  of  the  square  roughness  cases  and  the 
diamond  element  case  in  the  combined  pressure  gradient  region.  The  non-canonical  flows 
include  the  zero  and  favorable  pressure  gradient  diamond  element  cases. 

4.1  Canonical  Flows 

The  smooth  wall  and  square  roughness  Reynolds  shear  stress  profiles  are  compared  in  Fig. 
12a,  where  -  -uv  /  .  The  data  were  averaged  in  the  axial  direction  over  a  distance 

of  ~\.6-Sref.  The  vertical  lines  represent  the  scatter  in  the  shear  stress  data  over  the  averaged 
region;  as  indicated,  the  scatter  follows  the  uncertainty  estimates  in  Section  2.9.  The  smooth  and 
square  roughness  Reynolds  stresses  followed  the  expected  trends  from  the  literature  as  shown  in 
Fig.  lb;  recall  the  wall  shear  stress  was  estimated  from  the  turbulent  shear  stresses  via  the 
correlation  in  Fig.  lb.  The  corresponding  mean  deformations  across  the  boundary  layers  were 
computed  from  the  PIV  velocity  data.  The  deformations  were  normalized  by  the  root-mean- 
square  of  the  strain-rates  (^.  The  Bradshaw  (1974)  distortion,  individual  velocity  gradients, 
mean  dilatation  ( ),  strain-rate  and  spanwise  vorticity  were  tabulated  in  the  first  nine  columns 
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of  Table  4  for  y/S=  0.15,  0.4  and  0.6.  Also,  listed  in  column  10  are  the  Reynolds  shear  stresses  at 
the  same  locations.  The  remaining  columns  in  Table  4  are  discussed  in  Section  5. 

The  overall  structure  of  the  deformation  field  for  both  the  smooth  and  square  roughness  zero 
pressure  gradient  models  were  similar,  where  the  principal  strain  rate  dominated  and  the  extra 
strain-rates  and  dilatation  were  negligibly  small;  the  values  were  within  the  experimental  scatter. 
The  only  exception  was  a  moderate  near  wall  distortion  for  the  square  roughness  model  due  to 
the  presence  of  an  enhanced  dv  /  dy  gradient.  This  was  most  likely  the  result  of  the  relatively 
weak  waves  visible  in  Fig.  7a.  Hence,  the  observed  amplifications  of  the  turbulence  levels  were 
attributed  to  the  increased  form  drag  on  the  roughness  elements,  as  described  in  Goddard  (1959). 
Collectively,  the  small  dilatation,  weak  shocks  (Fig.  7b)  and  Reynolds  stress  scaling  in  Fig.  lb 
suggest  the  validity  of  Morkovin’s  hypotheses  for  the  square  roughness  model.  Moreover,  the 
shear  stress  and  deformation  results  suggest  that  J-type  roughness  effects  are  amenable  to 
numerical  simulation  through  modification  of  the  wall  boundary  condition  with  an  effective 
sand-grain  roughness  height;  e.g.,  the  k-comodc\  as  described  in  Wilcox  (2000). 
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Fig.  12  Reynolds  shear  stress  profiles  (/  Pre/U^ef  )i®  wind  tunnel  axes 


The  smooth  wall  turbulent  shear  stress  data,  with  the  imposed  curvature  induced  pressure 
gradients,  followed  previously  established  trends  from  the  literature  (Section  1.2),  where  the 
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favorable  pressure  gradients  had  the  expected  stabilizing  effect  on  turbulence  stresses.  As 
indicated  in  Fig.  12b,  the  weak  favorable  pressure  gradient  effects  were  relatively  modest,  where, 
in  the  lower  third  of  the  boundary  layer,  the  shear  stresses  were  attenuated  by  20%  with  respect 
to  the  zero  pressure  gradient  flow.  Relative  to  the  zero  pressure  gradient  values,  the  strong 
favorable  pressure  gradient  turbulent  shear  stress  values  were  reduced  by  about  40%  in  the  lower 
half  of  the  boundary  layer.  The  corresponding  strain-rates  are  listed  in  Table  4.  In  keeping  with 
Bradshaw’s  (1974)  definition  for  distortion,  the  values  for  d  are  presented  in  body  normal 
coordinates,  which  required  coordinate  system  rotations  of  -2.0°,  -4.1°  and  -9.7°  from  the  tunnel 
coordinate  system  for  the  weak  favorable  pressure  gradient  model  at  stations  2  and  3,  and  the 
strong  favorable  pressure  gradient  model  at  station  2,  respectively.  Imposing  the  weak  favorable 
pressure  gradient  resulted  in  a  0.1  distortion  strength  across  the  boundary  layer;  this  value  was 
considered  strong  based  on  Bradshaw’s  (1974)  criteria.  The <7v  /  cfy  velocity  gradient  was  the 
major  contributor  to  the  flow  distortion  across  the  boundary  layer.  For  the  smooth  wall  strong 
favorable  pressure  gradient  model,  the  total  distortion  values  in  the  outer-region  were  nominally 
0.3  -  0.4.  Flow  stretching  from  the  dilatation  remained  the  primary  contributor  to  as  dv  !  dy 
remained  dominant.  However,  the  du  /  dx  velocity  gradient  also  made  a  significant  contribution. 
The  dv !  dx  deformation  was  also  significant,  but  its  negative  value  resulted  in  a  25%  to  40% 
reduction  in  the  overall  flow  distortion.  Within  the  inner-region,  the  distortion  decreased  by  a 
factor  of  5  relative  to  outer-region  values,  principally  due  to  a  significant  decrease  in  dv !  dy . 
Luker  et  al.  (2000)  suggested  that  a  strengthened  dv !  dx  was  the  main  contributor  of  extra 
production  but  the  current  data  suggests  a  complimentary  effect,  where  all  of  extra  deformations 
tended  to  be  stabilizing. 
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Table  4.  Turbulence  Properties  across  the  Boundary  Layers 


y 

du 

du 

dv 

dv 

2k^ 

/  / 

u  u 

/  / 

VV 

)/l 

(PJ, 

iPxy)x 

s 

d 

dx 

dy 

dx 

dy 

^k,k 

^  xy 

pul, 

Ulf 

2k^ 

2 

2k^ 

2 

pul, 

2P 

2P 

2P 

Zero  Pressure  Gradient  -  Smooth 

0.15 

-0.01 

0.00 

1.00 

0.00 

-0.01 

-0.01 

0.50 

-0.50 

0.0009 

0.0066 

0.87 

0.13 

-10.0 

-0.4 

0.00036 

1.17 

-0.17 

-0.40 

0.40 

0.00 

-0.01 

1.00 

0.00 

0.00 

-0.01 

0.50 

-0.50 

0.0005 

0.0028 

0.88 

0.12 

-14.3 

-0.2 

0.00037 

1.16 

-0.16 

-0.15 

0.60 

0.02 

-0.01 

1.00 

0.00 

0.03 

0.02 

0.50 

-0.50 

0.0003 

0.0016 

0.88 

0.12 

-n.i 

-0.2 

0.00015 

1.17 

-0.17 

-0.20 

Zero  Pressure  Gradient  -  Square 

0.15 

0.04 

0.00 

1.00 

0.01 

0.03 

0.03 

0.51 

-0.50 

0.0020 

0.0382 

0.82 

0.18 

-4.1 

-0.4 

0.00259 

1.50 

-0.50 

-2.40 

0.40 

-0.02 

-0.02 

1.00 

0.00 

0.00 

-0.02 

0.50 

-0.50 

0.0016 

0.0072 

0.87 

0.13 

-12.7 

0.3 

0.00138 

1.16 

-0.16 

-0.24 

0.60 

-0.01 

-0.02 

1.00 

0.00 

0.00 

-0.02 

0.50 

-0.50 

0.0007 

0.0020 

0.85 

0.15 

-11.9 

0.3 

0.00019 

1.19 

-0.19 

-0.33 

Weak  Favorable  Pressure  Gradient 

-  Smooth 

0.15 

0.07 

0.01 

1.00 

0.00 

0.06 

0.07 

0.50 

-0.50 

0.0007 

0.0050 

0.87 

0.13 

-10.0 

-1.6 

0.00037 

1.29 

-0.29 

-0.62 

0.40 

0.07 

0.01 

1.00 

-0.01 

0.06 

0.07 

0.50 

-0.50 

0.0005 

0.0028 

0.86 

0.14 

-15.2 

-1.3 

0.00020 

1.29 

-0.29 

-0.27 

0.60 

0.10 

0.01 

0.99 

-0.04 

0.13 

0.14 

0.48 

-0.52 

0.0003 

0.0019 

0.85 

0.15 

-14.6 

-1.0 

0.00009 

1.45 

-0.45 

-0.21 

Weak  Favorable  Pressure  Gradient  -  Square 

0.15 

0.01 

0.00 

1.00 

0.00 

0.02 

0.02 

0.50 

-0.50 

0.0025 

0.0405 

0.92 

0.08 

-A.l 

-2.0 

0.00370 

1.12 

-0.12 

-0.83 

0.40 

0.04 

0.00 

1.00 

0.00 

0.04 

0.04 

0.50 

-0.50 

0.0016 

0.0088 

0.84 

0.16 

-17.6 

-1.6 

0.00115 

1.29 

-0.29 

-0.25 

0.60 

0.06 

0.00 

1.00 

-0.01 

0.08 

0.08 

0.50 

-0.50 

0.0007 

0.0049 

0.84 

0.16 

-16.0 

-1.3 

0.00044 

1.35 

-0.35 

-0.26 

Strong  Favorable  Pressure  Gradient  -  Smooth 

0.15 

0.08 

0.06 

0.99 

-0.06 

0.08 

0.14 

0.47 

-0.53 

0.0004 

0.0057 

0.74 

0.26 

-10.1 

-8.9 

-0.00018 

0.65 

0.35 

1.70 

0.40 

0.32 

0.05 

0.92 

-0.12 

0.36 

0.41 

0.40 

-0.52 

0.0003 

0.0023 

0.76 

0.24 

-14.8 

-7.6 

-0.00007 

-0.19 

1.19 

0.56 

0.60 

0.44 

0.04 

0.84 

-0.18 

0.51 

0.55 

0.33 

-0.51 

0.0003 

0.0017 

0.74 

0.26 

-12.4 

-6.3 

-0.00008 

0.06 

0.94 

0.21 

Strong  Favorable  Pressure  Gradient  -  Square 

0.15 

0.09 

-0.01 

0.99 

-0.04 

0.14 

0.13 

0.48 

-0.52 

0.0014 

0.0162 

0.74 

0.26 

-10.7 

-7.3 

-0.00025 

-0.80 

1.80 

5.12 

0.40 

0.17 

0.00 

0.97 

-0.06 

0.24 

0.24 

0.46 

-0.52 

0.0010 

0.0053 

0.77 

0.23 

-18.2 

-6.2 

0.00005 

5.84 

-4.84 

-3.43 

0.60 

0.26 

-0.01 

0.93 

-0.11 

0.36 

0.37 

0.41 

-0.52 

0.0006 

0.0028 

0.78 

0.22 

-13.8 

-5.0 

-0.00005 

-1.17 

2.17 

1.05 

Combined  Pressure  Gradient  -  Smooth 

0.15 

0.01 

-0.01 

1.00 

0.00 

0.02 

0.01 

0.50 

-0.50 

0.0015 

0.0091 

0.86 

0.14 

-12J 

-0.3 

0.00172 

1.21 

-0.21 

-0.30 

0.40 

-0.05 

-0.04 

1.00 

0.02 

-0.03 

-0.07 

0.51 

-0.49 

0.0005 

0.0018 

0.79 

0.21 

-27.9 

0.2 

0.00024 

1.29 

-0.29 

0.00 

0.60 

-0.29 

-0.06 

0.94 

0.10 

-0.31 

-0.37 

0.52 

-0.42 

0.0003 

0.0008 

0.73 

0.27 

-24.9 

-0.3 

0.00010 

1.10 

-0.10 

-0.35 

Combined  Pressure  Gradient  -  1 

Square 

0.15 

0.03 

0.03 

1.00 

0.00 

0.01 

0.04 

0.50 

-0.50 

0.0045 

0.0341 

0.82 

0.18 

-12.2 

-2.2 

0.00511 

1.36 

-0.36 

-0.66 

0.40 

0.02 

-0.01 

1.00 

0.01 

0.02 

0.01 

0.51 

-0.50 

0.0021 

0.0073 

0.83 

0.17 

-21.7 

-1.2 

0.00130 

1.26 

-0.26 

-0.07 

0.60 

-0.12 

-0.02 

0.99 

0.06 

-0.15 

-0.17 

0.53 

-0.47 

0.0010 

0.0025 

0.74 

0.26 

-33.5 

-1.4 

0.00041 

1.30 

-0.30 

-0.06 

Combined  Pressure  Gradient  -  Diamond 

0.15 

-0.03 

-0.05 

1.00 

0.00 

0.02 

-0.03 

0.50 

-0.50 

0.0025 

0.0389 

0.91 

0.09 

-4.4 

-0.6 

0.00852 

1.08 

-0.08 

-0.37 

0.40 

-0.02 

-0.04 

1.00 

0.03 

-0.01 

-0.05 

0.52 

-0.49 

0.0007 

0.0048 

0.82 

0.18 

-21.5 

-0.2 

0.00095 

1.23 

-0.23 

-0.10 

0.60 

-0.30 

-0.06 

0.94 

0.09 

-0.31 

-0.37 

0.52 

-0.43 

0.0003 

0.0015 

0.75 

0.25 

-29.9 

-1.2 

0.00030 

1.13 

-0.13 

-0.20 

Notes:  (1)  Columns  du  /  dx  through  are  normalized  by  g  ;  (2)  the  subscript  X  indicates  in  principal  axes. 
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The  strong  favorable  pressure  gradient  veloeity  eases,  listed  in  Table  4,  were  the  only  ones 
that  had  a  significant  impact  on  the  bulk  dilatation  and  the  Sxy  strain-rate.  For  all  of  the  other 
cases,  the  dilation  was  negligible,  and  Sxy  was  nominally  0.5.  For  the  smooth  strong  favorable 
pressure  gradient,  the  trend  for  the  bulk  dilatation  was  that  the  magnitude  became  significantly 
larger  with  pressure  gradient,  and  for  the  strong  favorable  pressure  gradient,  the  mean  dilatation 
was  approximately  30  -  50%  of  the  principal  strain  rate.  The  Sxy  strain-rate  decreased  from  0.47 
to  0.33  moving  from  the  inner  to  outer  regions  of  the  boundary  layer.  From  an  eddy  viscosity 
modeling  perspective,  this  would  indicate  a  decrease  in  the  Reynolds  shear  stress,  with  the  most 
pronounced  effect  in  the  outer  layer.  When  the  present  strong  favorable  pressure  gradient  data  in 
Fig.  12b  were  expressed  in  body  normal  coordinates  as  described  in  Luker  et  al  (2000),  the  shear 
stress  levels  in  the  outer  region  were  significantly  reduced,  which  is  consistent  with  the  mixing 
length  model.  It  is  expected  that  the  pressure  work  terms,  which  do  not  have  incompressible 
counterparts,  contributed  to  the  turbulent  shear  stress  transport  dynamics.  More  specifically,  the 
ratio  of  the  pressure-work  to  production  for  the  Reynolds  shear  stress  was  crudely  estimated 
using  the  strong  Reynolds  analogy  [Smits  and  Dussauge  (1996)].  Accounting  for  axial  pressure 
gradients  and  the  primary  production  component  resulted  in  the  following  expression; 

P„/  P - [(7- V}M^]  /  Tyy)  {pJJ^  / pU)  [{p  -p^l \ pJJ]A  {S I L).  For  the  present  strong 

favorable  pressure  gradient  model,  the  first  term  in  this  expression  was  ~4.0,  the  second  two 
were  both  order  1.0,  the  pressure  coefficient  was  0.07  and  the  last  term  was  ~0.5.  Thus,  the 
pressure- work  to  production  ratio  was  roughly  a  modest  15%. 

Coupling  the  square  roughness  with  the  favorable  pressure  gradients  introduced  a 
competition  between  the  destabilizing  roughness  and  stabilizing  pressure  gradients.  For  example, 
the  smooth  wall  and  square  surface  roughness  model  distortions  in  the  strong  favorable  pressure 


36 


gradient  region  were  similar  in  the  inner-boundary  layer  region.  However,  the  eontributions  from 
the  individual  gradients  differed.  Specifically,  the  rough  wall  axial  deformations,  including  the 
stabilizing  dv !  dx  were  significantly  reduced,  while  dv !  dy  increased  by  75%.  In  the  outer- 
region,  the  square  roughness  distortion  was  40%  lower  than  the  corresponding  smooth  wall 
value.  The  net  effect  was  reduced  changes  in  the  mean  dilatation  and  Sxy.  The  results  of  these 
competing  effects  are  apparent  by  comparing  the  square  roughness  Reynolds  shear  stress  surveys 
in  Figs.  12a  and  12b.  Near  the  wall,  the  strong  favorable  pressure  gradient  Reynolds  stress  values 
were  similar  to  the  zero  pressure  gradient  levels.  However,  further  out  in  the  boundary  layer,  a 
decline  in  turbulent  shear  stress  relative  to  the  zero  pressure  gradient  model  was  observed,  with 
reductions  on  the  order  of  50%.  These  combined  effects  require  a  reasonable  degree  of  inherent 
fidelity  in  a  turbulence  model  to  duplicate  the  intricate  Reynolds  stress  strain  rate  relationship  for 
this  case.  Moreover,  the  pressure-work  to  production  ratio  is  expect  to  have  been  enhanced  by 
the  roughness  as  Latin  and  Bowersox  (2000)  have  shown  that  the  correlation  between  the  density 
and  velocity  fluctuations  increased  with  increasing  roughness  height.  Thus,  the  pressure-work 
was  considered  important  in  generating  turbulence  for  this  flow.  However,  the  exact  role  was  not 
considered  any  further  as  direct  measurements  were  not  achieved  in  this  study.  Overall,  even 
with  these  shortcomings,  it  is  believed  that  current  two-equation  turbulence  models  that  treat  this 
roughness  topology  via  a  modified  boundary  condition  will  capture  the  overall  trends  associated 
with  the  competing  effects  provided  that  appropriate  compressibility  corrections  are  included. 

The  shear  stress  profiles  from  the  combined  pressure  gradient  region  are  given  in  Fig.  12c. 
The  changes  in  the  turbulent  stress  profile  for  the  combined  pressure  gradient  region  indicated 
the  turbulent  stress  was  dependent  on  the  proximity  to  the  dual  sources  of  flow  strain;  the 
principal  strain-rate  and  the  recompression.  The  corresponding  velocity  gradient  data  are  listed  in 
Table  4.  The  sign  of  all  extra  strain-rates  were  reversed  as  compared  to  the  favorable  pressure 
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gradient  location.  This  resulted  in  a  negative  overall  distortion.  The  distortion  values  near  the 
wall  were  much  weaker  than  those  in  the  outer  region,  where  in  the  outer  region,  flow  turning 
from  the  recompression  shock  which  was  more  pronounced.  In  the  inner  region,  the  bulk 
dilatation/compression  was  essentially  zero  small.  However,  in  the  outer  region,  the  flow 
elements  were  significantly  compressed  and  the  vorticity  was  reduced. 

The  Reynolds  shear  stresses  for  the  rough  wall  combined  pressure  locations  are  also  shown 
in  Fig.  12c.  The  square  roughness  and  adverse  pressure  gradient  coupled  to  produce  turbulent 
shear  stress  levels  that  were  400%  above  the  smooth  wall  values.  The  profile  shape  for  this  case 
also  had  a  peak  that  occurred  near  yl5=  0.3.  For  the  combined  pressure  gradient  diamond 
roughness  model,  the  profile  shape  for  the  turbulent  shear  stress  was  similar  to  that  of  the  smooth 
and  square  roughness  case;  however,  the  levels  were  nominally  40%  lower  than  the  square 
roughness  model.  These  lower  turbulence  levels  were  expected,  as  the  skin  friction  at  the  wall 
was  primarily  a  function  of  the  roughness  element  drag  [Liepmann  and  Goddard  (1957)]  and  the 
drag  coefficient  for  the  streamlined  diamond  topology  was  expected  to  be  lower  than  that  for  the 
blunt  square  shapes.  Also  implicit  in  this  argument  are  the  assumptions  that  Townsend’s  scaling 
(1976)  holds  and  the  turbulence  levels  responded  to  the  local  surface  condition.  Adding 
roughness  had  minimal  effects  on  the  compression,  Sxy  and  vorticity  for  the  combined  model  as  it 
did  for  the  favorable  pressure  gradient  models. 

4.2  Non-Canonical  Flows 

The  diamond  element  surface  roughness  models  for  the  zero  and  favorable  pressure  gradient 
cases  produced  Reynolds  shear  fields  that  had  a  systematic  dependence  on  the  relative  location  to 
the  surface  topology.  To  demonstrate  the  nature  of  the  flow,  contours  of  the  Reynolds  shear 
stress  and  the  dv/dy  gradient  for  the  zero  and  strong  favorable  pressure  models  presented  in  Fig. 
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13.  To  better  quantify  the  flow  behavior,  axial  plots  of  the  deformation  and  distortions  aty/S= 
0.15,  0.4  and  0.6  are  given  in  Fig.  14;  the  symbol  legend  is  in  the  right-most  plot  in  Fig.  14a.  The 
left  ordinate  scale  corresponds  to  the  distortion,  and  the  strain-rates,  normalized  by  the  total 
strain  rate,  are  on  the  right.  The  weak  favorable  pressure  gradient  cases  were  excluded  from  the 
present  discussion,  as  their  results  can  be  inferred  from  the  trends  presented  in  Table  4  and  Fig. 

14.  For  the  zero  pressure  gradient  model  aty/S=  0.15,  the  deformations  appear  erratic.  This  was 
not  surprising  as  the  flow  near  the  leading  edges  of  the  roughness  elements  showed  signs  of  a 
localized  disturbance  (see  Fig.  11).  Additionally,  laser  wall  reflections  had  a  greater  influence 
closer  to  the  wall  and  possibly  affected  the  results  (See  Fig.  6). 

Traversing  from  left-to-right,  for  the  zero  pressure  gradient  flow  in  Fig.  13,  the  flow  crossed 
the  roughness  element  leading-edge  shocks,  which  caused  a  band  of  negative  dvjdy  and 
increased  turbulence  shear  stresses.  The  leading-edge  shear  stress  increases  were  the  result  of 
shock  turbulence  interactions  and  the  subsequent  adverse  pressure  gradient  on  the  front  half  of 
the  roughness  element  (see  Fig.  8).  Further  downstream,  the  strain  rate  changed  sign  and  the 
turbulence  levels  decreased.  This  decrease  was  the  result  of  the  local  favorable  pressure  gradients 
on  the  back  half  of  the  roughness  elements  (see  Fig.  8).  Therefore,  in  the  axial  direction,  the 
Reynolds  shear  stress  had  a  periodic  variation,  where  the  minima  values  were  on  the  order  of  the 
smooth  wall  models,  while  maxima  values  were  doubled.  The  rapid  response  of  the  Reynolds 
shear  stress  to  the  strain-rate  field  indicated  that  linear  effects  were  predominant  in  the  transport 
dynamics.  The  corresponding  distortions,  in  Fig.  14,  varied  between  ±0.4,  with  the  negative 
distortions  corresponding  to  the  adverse  pressure  gradient  just  after  the  roughness  element 
leading  edge,  while  the  positive  distortions  resulted  from  the  favorable  pressure  gradient 
generated  by  the  roughness  element  trailing  edge  expansions.  The  distortion  behavior  was 
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strongly  coupled  to  dv  !  dy  ,dLS,  du  I  dx  and  dv  !  dx  tended  to  cancel.  The  mean  dilatation  and  Sxy 
tracked  the  distortion,  and  the  vorticity  was  nominally  constant  at  -0.5. 
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Fig.  13  Strain-rate  and  Reynolds  shear  stress  fields  for  the  zero  (left  images)  and  strong  favorable  (right 
images)  pressure  gradient  models  with  diamond  roughness. 


For  the  strong  favorable  pressure  gradient  model,  the  global  pressure  gradient  coupled  with 
the  local  gradients  generated  by  diamond  roughness  elements  to  produce  a  highly  stratified  (in 
the  axial  direction)  shear  stress  field;  see  Fig.  13.  As  with  the  zero  pressure  gradient  case. 
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relatively  large  shear  stress  levels  were  observed  near  the  leading  edge  shoek.  In  between  the 
shock,  the  turbulence  was  highly  damped  by  the  presence  of  the  global  and  local  expansions, 
where  the  Reynolds  stresses  were  reduced  to  values  below  those  of  the  zero  pressure  gradient 
smooth  wall  boundary  layer.  The  large  shear  stress  levels  at  the  leading  edge  demonstrated  that 
shock  turbulence  interaction  was  a  dominant  shear  stress  producing  feature  in  this  flow.  The  data 
in  Fig.  14b  demonstrated  that  the  global  favorable  pressure  gradient  amplified  all  distortions, 
both  positive  and  negative.  For  example,  the  distortions  varied  from  -0.75  to  1.45.  The  mean 
dilatation  peaks  varied  from  -0.5  to  1.0.  The  Sxy  strain  rate  varied  from  0.25  to  0.45,  where  the 
minimum  corresponded  to  the  peak  distortion,  and  the  vorticity  was  nominally  constant  at  -0.45. 
These  changes  likely  were  primarily  due  to  a  decrease  in  the  overall  strength  of  the  principal 
strain-rate  {du  !  dy ),  and  enhanced  extra  strain-rates  due  to  the  orientation  of  the  wall  normal 
with  the  shock  and  expansion  waves.  Supporting  evidence  comes  from  the  schlieren  image  in 
Fig.  7b,  where  the  shock  waves  generated  by  the  diamond  element  surface  roughness  became 
almost  normal  to  the  wall,  and  thus  contributed  little  to  the  wall  normal  principal  strain-rate 
while  enhancing  the  extra  velocity  gradients. 

The  data  in  Fig.  13  indicated  that  the  Reynolds  shear  stress  responded  strongly  to  the  local 
distortions,  which  implied  that  the  local  effects  appeared  to  prevail  over  history  effects.  The 
corresponding  intricate  Reynolds-stress/strain-rate  relationships  imply  that  numerical  simulations 
must  resolve  the  flow  structure  within  the  roughness  topology  to  accurately  capture  the  flow 
structure.  The  pressure  work  terms  were  expected  to  be  more  important  for  this  case,  as 
compared  to  the  global  gradient  cases,  because  similar  pressure  changes  were  observed  over  a 
smaller  spatial  extent;  6IL  ~  1.  As  was  the  case  in  the  previous  section,  the  exact  role  was  not 
considered  any  further  as  direct  measurements  were  not  achieved.  However,  it  was  observed  that 
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magnitude  of  the  local  Reynolds  stress  response  to  the  local  strain  rates  followed  the  pattern 
indicated  by  the  canonical  flows.  This  result  suggests  that  current  two-equation  turbulence 
models  with  suitable  compressibility  corrections  should  capture  the  local  trends. 

4.3  Section  Summary 

The  Reynolds  stress  response  to  the  canonical  flow  deformations  behaved  in  an  expected 
manner.  For  the  zero  pressure  gradient  square  roughness  case,  compressibility  effects  were 
minimal  and  the  flow  structure  across  the  boundary  layer  appeared  to  scale  with  the  roughness 
surface  condition  as  described  by  Townsend  (1976).  Hence,  this  case  appears  to  be  amenable  to 
current  prediction  methods.  The  Reynolds  stress  response  to  the  local  strain  rates  followed  the 
expected  trend,  where,  as  indicted  by  Bradshaw  (1974),  the  distortions  affected  the  flow  by  an 
order  of  magnitude  more  than  would  be  anticipated  by  the  associated  extra  production.  This 
result  suggests  that  current  two-equation,  or  higher,  turbulence  models  with  suitable 
compressibility  corrections  may  capture  the  correct  trends.  Collectively,  the  results  in  Section  4. 1 
suggest  that  modem  models  may  capture  the  competing  roughness  and  pressure  gradient  effects 
for  the  canonical  cases  with  square  roughness. 

The  diamond  element  roughness  produced  local  flow  phenomena  that  altered  the  local  state 
of  the  boundary  layer,  representing  a  marked  departure  from  the  canonical  flows.  The  intricate 
Reynolds-stress/strain-rate  relationship  implies  that  numerical  simulations  should  resolve  the 
flow  within  the  roughness  topology  to  accurately  capture  the  overall  stracture.  The  magnitude  of 
the  local  Reynolds  stress  response  to  the  local  strain  rates  followed  the  pattern  indicated  by  the 
canonical  flows.  This  result  suggests  that  current  turbulence  models  with  suitable  compressibility 
corrections  may  capture  the  local  trends.  These  findings  motivate  the  examination  of  the 
turbulence  production  and  evolution  mechanisms  discussed  in  Section  5. 
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(b)  Strong  Favorable  Pressure  Gradient 
Fig.  14  Axial  profiles  of  diamond  roughness  model  flow  distortions 
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5.  Discussion  of  the  Turbulence  Evolution  and  Production 


The  intricate  Reynolds  stress/strain-rate  relationships  described  in  the  previous  section  for 
both  the  canonical  and  non-canonical  flows  motivated  a  closer  examination  of  the  turbulence 
production  and  evolution  mechanisms.  To  help  facilitate  the  discussion,  an  idealized 
representation  of  turbulence  stress  orientation  evolution  for  the  strong  pressure  gradient  model  is 
shown  in  Fig.  15.  The  ordinate  for  the  wall  coordinate  system,  is  represented  by  the  arrow 
originating  along  the  floor,  which  for  the  F*  and  3'^'*  measurement  locations  correspond  to  the  test 
section  coordinate  system.  The  dashed  axes  at  each  location  represent  the  principal  axes  for  the 
turbulent  stress  aty/S=  0.4.  The  turbulent  stress  turning  angle  a  is  the  rotation  angle  from  the 
test  section  coordinate  system  to  the  turbulence  stress  principal  axes.  The  ellipsoids,  illustrate  the 
allocation  and  orientation  of  the  principal  stresses.  The  test  section  coordinate  system  is  shown 
for  reference  at  the  2"^*  measurement  location.  Similarly  at  the  second  test  location,  the  flow- 
turning  angle  (j)  represents  the  rotation  angle  from  the  test  section  coordinate  system  to  the  axis 
along  the  mean  velocity  vector.  This  was  omitted  from  the  axes  at  locations  1  and  3  since,  as  will 
be  seen  later,  the  mean  velocity  vector  was  nearly  aligned  with  the  axial  direction. 

The  turbulence  stress  principal  axes  were  not  known  a  priori.  Thus,  the  data  were  collected  in 
the  test  section  coordinate  system  and  transformed  to  their  principal  stresses  via  the  coordinate 

transformation  given  in  equations  5.1a-c,  with  the  condition  that  'v^ '  is  equal  to  zero. 

'u^ '  =  w'u'cos^  C(r-i-2i/'v'cosc(rsina-i-v'v'sin^  a  (5.1a) 

'v^'  =  u'w'sin^  ctr-2w'v'cosc(rsina-i-v'v'cos^  a  (5.1b) 

~  (w'w'-v'v')cosc(rsina-i-i/'v'(cos^  a-i-sin^  a)  (5.1c) 
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For  the  present  study,  the  measured  Ic  values  were  defined  as  ^{u'u'+v'v') .  Neglecting  the 
lateral  velocity  gradients,  the  production  terms  were  calculated  in  principal  axes  as  follows: 

(5.2a) 

(/>„),  =2(T„),(3v/Sti.),  (5.2b) 

(PJz  ==  (5-2c) 

The  production  of  turbulent  kinetic  energy  was  normalized  by  Idj pUl^  ,  while  the  production 
components  were  normalized  by  2P. 


Summarized  in  Table  4  (columns  11  -  19)  are  the  turbulence  kinetic  energy,  normal  stress 
components,  turbulent  stress  turning  angle  a,  flow  turning  angle  (p,  the  turbulent  kinetic  energy 
production  and  the  production  components  for  the  canonical  cases  at  three  transverse 
measurement  locations  {y/S=  0.15,  0.4  and  0.6).  As  the  weak  favorable  pressure  gradient  data 
were  only  slightly  altered  as  compared  to  the  zero  pressure  gradient  data,  they  are  omitted  from 
the  subsequent  discussion;  these  data  are  included  in  Table  4  for  completeness.  The 
corresponding  non-canonical  flow  results  are  presented  as  line  plots  parallel  to  the  roughness 
elements  at  the  same  transverse  locations  as  listed  in  Table  4. 
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5.1  Canonical  Flows 


As  indicated  by  comparison  of  the  first  and  second  rows  in  Table  4,  adding  the  square 
roughness  topology  led  to  signifieant  inereases  in  the  near  wall  turbulent  kinetie  energy,  where  a 
480%  inerease  was  observed  aty/<^=  0.15.  Increased  momentum  loss  at  the  roughness  elements 
due  to  form  drag  is  believed  to  be  the  major  driver  for  the  dramatie  increase  in  With 
increasing  distance  from  the  wall  however,  the  relative  increase  in  was  reduced.  Wall 
roughness  also  signifieantly  inereased  the  relative  amount  of  the  transverse  kinetie  energy  near 

the  wall.  Specifieally,  aty/^=  0.15,  18%  of  was  due  to  'v^ '  for  the  rough  wall,  where  the 

corresponding  smooth  wall  value  was  13%,  indicating  an  increase  in  secondary  turbulent  flow 
structures  near  the  roughness  elements.  The  turbulent  stress  turning  angle  (a)  for  the  zero 
pressure  gradient  smooth  wall  model  were  rotated  about  -12°  from  the  test  seetion  eoordinate 
system  while  the  loeal  flow  angle  (j)  remained  parallel  with  the  wall,  which  agrees  well  with  the 
smooth  wall  zero  pressure  gradient  LDV  data  from  Luker  et  al.  (2000)  and  Latin  and  Bowersox 
(2000).  While  the  orientation  in  the  outer  portion  of  the  boundary  layer  was  similar  for  both  wall 
surfaces,  a  rotated  significantly  from  -10°  to  -4.1°  for  the  rough  wall  case  aty/^=  0.15.  The 
redistribution  of  the  turbulenee  stress  tensor  away  from  shear  and  towards  eompression  helps  to 
explain  the  unequal  modifieation  of  the  turbulence  stress  measurements  in  the  wind  tunnel  axes 
system,  where  the  turbulence  shear  stress  doubled  and  the  kinetic  energy  increased  five-fold. 

The  overall  increase  in  produetion  (P)  for  the  square  roughness  zero  pressure  gradient  model 
was  in  line  with  observed  increases  in  k^  at  all  locations.  The  composition  of  P  from  {Pxx)x  and 
{Pyy)i  remained  constant  throughout  the  boundary  layer  for  the  smooth  wall  zero  pressure 
gradient  model,  with  produetion  due  primarily  to  the  {Pxx)x  term,  while  the  {Pyy)x  term  resulted  in 
a  slight  reduetion  in  overall  P.  The  eomposition  of  P  for  the  square  roughness  zero  pressure 
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gradient  model  at  the  outer  two  locations  was  similar.  However,  the  near  wall  {Pxx)i  was  larger 
and  the  stabilizing  effect  from  {Pyy)x  was  also  more  substantial.  The  square  roughness  inner- 
region  {Pxy)x  term  was  ~500%  greater  than  the  corresponding  smooth  wall  value,  which  indicated 
the  roughness  elements  resulted  in  a  highly  anisotropic  production  tensor.  Within  the  inner- 
region  {Pxy)x  had  a  magnitude  that  was  240%  greater  than  the  cumulative  contributions  from  the 
normal  production  components.  Thus,  in  this  region  the  production  tensor  principal  axes  were 
misaligned  with-respect-to  the  turbulence  stress  principal  axes.  These  results  are  consistent  with 
the  idea  that  the  principal  near  wall  production  mechanisms  are  altered  from  the  conventional 
bursting  process  typically  seen  for  smooth  wall  boundary  layers  towards  localized  production 
driven  by  roughness  element  induced  flow  structures  such  as  wake  vorticies  [Grass  (1971)].  The 
increase  in  {Pxy)x  near  the  wall  is  also  relevant  to  the  pressure-strain  redistribution,  where  for 
example.  Launder  et  al  (1975)  argued  the  essential  features  of  the  rapid  part  of  the  pressure-strain 
is  linearly  proportional  to  the  production  anisotropy  for  low-speed  flows.  The  smallness  of  the 
dilatation  for  this  case  suggests  applicability.  The  additional  data  listed  in  Table  4  may  provide 
additional  insight  into  the  rapid  pressure-strain  via  more  general  models.  However,  such  efforts 
were  not  attempted  here,  where,  instead,  the  focus  was  on  directly  measured  quantities. 
Nonetheless,  the  observed  modification  of  the  pressure-strain  redistribution  may  be  important  to 
accurately  predict  the  turbulence  stress  component  via  a  seven  equation  Reynolds  stress  transport 
model;  this  implication  applies  to  all  of  the  remaining  canonical  flow  cases  described  below. 

The  curvature  induced  strong  favorable  pressure  gradient  had  the  expected  effect  of  reducing 
the  overall  for  the  smooth  wall  cases,  where  the  percent  changes  were  -14%,  -18%  and  6.2% 
at  y/d=  0.15,  0.4  and  0.6,  respectively.  The  turbulence  within  the  strong  favorable  pressure 
gradient  boundary  layer  was  significantly  more  isotropic  than  the  smooth  zero  pressure  gradient 
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case,  where  25%  of  was  due  to  v^'v^'  and  75%  from  u^'u^' .  This  result  agrees  with  the 

intermittency  conclusions  from  Luker  et  al.  (2000),  where  it  was  also  concluded  the  turbulence 
length  scales  were  reduced,  and  thus  the  fluctuation  were  more  readily  dissipated.  The  axial 
normal  stress  component  production  was  significantly  reduced,  with  a  negative  value  near  y/S = 
0.4,  which  indicated  that  turbulence  energy  was  being  transferred  back  to  the  mean  flow  [Arnette 
et  al.  (1996)].  The  positive  sign  on  the  transverse  stress  production  explains  the  relative  increase 
in  transverse  shear  stress  as  compared  to  the  zero  pressure  gradient  flow.  In  the  strong  favorable 
pressure  gradient  flow,  the  positive  transverse  production  meant  that  this  component  was  being 
directly  energized  by  the  mean  flow,  where  for  the  zero  pressure  gradient  case,  the  negative 
transverse  production  results  in  an  energy  flow  from  the  turbulence  to  the  mean  flow,  and  the 
transverse  energy  fluctuations  are  the  result  of  the  cross  coupling  between  the  axial  and 
transverse  components  (e.g.,  pressure-strain  redistribution).  The  sign  change  of  the  shear  stress 
production  explains  the  overall  reduction  in  the  shear  stress  plots  in  Fig.  12b.  More  intuitively, 
this  shear  reduction  can  be  explained  by  again  referring  to  Fig.  15  and  the  values  of  a  in  Table  4. 
For  the  strong  favorable  pressure  gradient  model  the  difference  between  the  principal  stresses 
and  the  main  flow  direction  (i.e.,  a  -  (p)  was  -1.2°,  -7.2°  and  -6.1°  aiy/d=  0.15,  0.4  and  0.6 
respectively.  This  separation  was  much  less  significant  than  the  values  of  -9.6°,  -14.1°  and  - 
12.5°,  observed  for  the  zero  pressure  gradient  smooth  model.  Thus,  the  compressive  stresses 
exerted  on  the  flow  were  enhanced  while  the  shear  was  minimized.  This  was  the  direct  result  of 
the  extra  strain  rates  altering  the  production,  and  hence  structure  of  the  turbulence. 

As  with  the  zero  pressure  gradient  models,  adding  surface  roughness  to  the  strong  favorable 
pressure  gradient  models  significantly  increased  relative  to  the  corresponding  smooth  wall 
models.  Specifically,  comparing  the  strong  favorable  pressure  gradient  rough  and  strong 
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favorable  pressure  gradient  smooth  wall  results,  increases  in  of  185%,  130%  and  65%,  foxy/5 
=  0.15,  0.4  and  0.6  were  observed,  respectively.  This  result  suggests  a  competition  between  the 
global  pressure  gradient  and  the  flow  structure  established  by  the  roughness  topology. 
Specifically,  the  roughness  elements  produce  increased  shear  stress  through  form  drag,  and  the 
favorable  pressure  gradient  was  working  to  damp  the  turbulence  through  the  mechanisms 
described  earlier.  Hence,  even  though  the  strong  favorable  pressure  gradient  square  roughness 
turbulence  levels  increased  relative  to  the  smooth  strong  favorable  pressure  gradient  model,  the 
levels  were  significantly  lower  than  zero  pressure  gradient  square  roughness  cases.  The  overall 
production  and  production  components  demonstrate  this  competition.  Adding  the  roughness 
tended  to  slightly  decrease  the  stress  principal  axis  direction,  with  the  largest  change  of  -3.4° 
near  y/5=  0.4.  The  anisotropy  for  the  rough  strong  favorable  pressure  gradient  case  was  similar 
to  that  of  smooth  strong  favorable  pressure  gradient  case. 

In  the  combined  pressure  gradient  region  of  the  smooth  wall  model,  was  38%  larger  than 
the  corresponding  zero  pressure  gradient  values  at  y/5  =  0.15.  However,  for  the  outer  two 
locations,  the  turbulence  levels  had  not  yet  been  amplified  back  to  either  the  zero  pressure 
gradient  or  strong  favorable  pressure  gradient  levels.  This  is  understandable  as  the  production  of 
at  y/5  =  0.15  was  about  an  order  of  magnitude  larger  than  in  the  outer  regions.  Thus,  the 
recovery  of  was  slower  in  the  outer  regions.  The  lower  outer  region  levels,  relative  to  the 
upstream  strong  favorable  pressure  gradient  values,  imply  that  k^  continued  to  decrease  axially 
beyond  the  strong  favorable  pressure  gradient  measurement  location.  The  structural  content  of 
the  turbulence  at  y/5=  0.15  was  similar  to  that  for  the  zero  pressure  gradient  model.  However,  in 
the  outer  regions,  the  turbulence  appeared  to  be  more  isotropic  and  the  stress  tensor  principal 
axes  angles  were  significantly  lower.  These  trends  are  shown  schematically  in  Fig.  15,  where  it 
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can  be  seen  that  the  turbulence  went  from  highly  anisotropic  turbulence  towards  a  more  isotropic 
with  a  favorable  pressure  gradient,  and  back  towards  anisotropy  once  passing  through  the 
adverse  pressure  gradient.  The  production  component  values  were  similar  to  those  from  the  zero 
pressure  gradient  case,  and  hence  drove  the  turbulence  back  towards  zero  pressure  gradient 
values.  The  trends  for  the  two  rough  wall  combined  pressure  gradient  models  were  similar  to  the 
smooth  wall  results  with  the  addition  of  the  enhanced  production  due  to  the  roughness.  The 
higher  production  for  the  diamond  element  model  at  y/^=  0.15  was  attributed  to  the  thinner 
boundary  layers  and  steeper  velocity  gradients  in  this  region. 

5.2  Non-Canonical  Flows 

As  noted  previously  in  Section  4.1,  the  diamond  element  surface  roughness  created  strain 
rate  fields  that  locally  altered  turbulent  stress  phenomena.  Contours  of  the  Reynolds  shear  stress 
production,  in  test  section  coordinates,  for  the  zero  pressure  gradient  and  strong  favorable 
pressure  gradient  models  are  displayed  in  Fig.  16.  The  values  were  normalized  by  the  reference 
length  scale  divided  by  the  mean  density  and  the  cube  of  the  reference  velocity  /  pUl^f).  In 

both  images,  there  are  black  structures  near  the  lower  surface.  These  structures  were  the  result  of 
selecting  the  grayscale  range  to  highlight  the  changes  in  the  outer  region  boundary  layer  (y/S> 
0.2).  The  near  wall  production  Pxy  was  comparatively  much  larger  due  to  the  strong  principal 
strain-rate  and  elevated  turbulence  levels  due  to  the  roughness  elements.  Both  the  zero  and 
strong  favorable  pressure  gradient  model  showed  similar  structure  in  the  Pxy  contours.  The 
systematic  positive  and  negative  production  regions  are  strikingly  periodic.  The  dark  regions 
(contour  levels  ~  -0.002)  corresponded  to  increased  production  of  the  Reynolds  shear  stress  due 
to  the  leading  edge  shocks  and  the  subsequent  compressive  strain-rate  field  over  the  front  half  of 
the  roughness  element.  In  the  light  regions  between  the  shocks,  the  production  changed  sign  as 
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the  flow  experienced  stabilizing  favorable  pressure  gradients  along  the  aft  half  of  the  roughness 
element.  The  strong  favorable  pressure  gradient  model  produced  stronger  shocks  and  larger 
regions  of  stabilization  between  the  shocks.  This  result  confirms  that  the  shock  turbulence 
interactions  were  responsible  for  part  of  the  turbulence  level  increases  near  the  leading  edge  of 
the  roughness  element  where,  as  described  in  Mahesh  et  al.  (1997),  the  turbulence  amplification 
across  the  shock  depends  on  the  angle  between  the  shock  and  the  shear  layer.  Hence,  the 
amplifications  varied  across  the  boundary  layer  (see  the  curved  shocks  in  Fig.  7b)  and  axially  for 
the  pressure  gradient  cases. 


(a)  ZPG  Model  (b)  SFPG  Model 

Fig.  16  Diamond  roughness  element  model  turbulent  shear  stress  production 


Line  plots  of  the  turbulence  structure  are  given  in  Fig.  17  to  better  quantify  the  periodic 
nature  of  flow  structure.  The  legend  is  in  the  right-most  plot  in  Fig.  17a.  The  weak  favorable 
pressure  gradient  cases  were  excluded,  as  the  trends  were  similar  to  zero  pressure  gradient  case. 
For  the  zero  pressure  gradient  data  in  Fig.  17a,  and  a  was  observed  to  vary  significantly  aty/S 
=  0.15.  These  variations  were  due  to  the  flow  disturbances  previously  mentioned  near  the 
element  leading  edges.  Although,  the  y/S  =  0.15  data  are  included  for  completeness,  the 
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subsequent  discussion  will  be  limited  to  the  outer  two  locations.  Focusing  then  on  the  data  in  Fig. 
17a  aty/S=  0.4  and  0.6,  the  overall  average  value  of  was  close  to  the  smooth  wall  values,  with 
a  variation  of  ±20%  and  ±10%  at  j/<?=  0.4  and  0.6,  respectively.  The  sign  of  the  variation 
depended  on  whether  the  datum  point  was  in  the  local  adverse  or  favorable  pressure  gradient 
region  of  the  flow  established  by  the  surface  element.  However,  the  allocation  of  turbulent  stress 

was  nearly  constant  in  the  axial  direction,  with  85%  of  contained  within  for  both 

locations.  The  principal  stress  axis  angle  exhibited  strong  variations,  between  -10°  to  -20°  in  the 
favorable  and  adverse  pressure  gradient  regions,  respectively.  On  the  present  scale,  the  flow- 
turning  angle  appeared  to  be  essentially  constant.  However,  referring  back  v  717^/ contour  (Fig. 
10),  the  flow  turning  angles  were  systematic  as  described  in  Section  3.5. 

Axial  profiles  of  the  production  for  the  zero  pressure  gradient  and  strong  favorable  pressure 
gradient  models  are  presented  in  Fig.  18  for  y/S=  0.15,  0.4  and  0.6.  The  nondimensionalization 
follows  that  in  Table  4.  For  each  plot,  the  production  of  turbulent  kinetic  energy  is  on  the  left 
ordinate,  while  the  production  of  the  individual  stresses  is  on  the  right.  The  legend  is  in  the 
rightmost  plot  for  the  zero  pressure  gradient  case.  In  the  compression  region  for  y/S  =  0.4  and 
0.6,  the  scalar  production  P  had  a  maximum  value  that  was  similar  in  magnitude  to  that  observed 
for  the  square  surface  roughness  model  at  the  same  boundary  layer  locations  (in  Table  4).  In  the 
expansion  regions,  small  pockets  of  negative  production  similarly  existed.  The  (Pxx)^.  production 
term  was  dominant  in  the  leading  edge  compression  regions.  Conversely,  the  (Pyy)x  term  was 
dominant  in  the  areas  of  negative  production.  Finally,  the  (Pxy)x  term  remained  nearly  constant  at 
about  40%  of  P  at  all  boundary  layer  locations.  The  observed  production  fields  indicated  an 
involved  balance  between  production  and  redistribution,  and  this  balance  responded  rapidly  to 
the  local  flow  distortions.  It  is  likely  that  these  intricate  processes  will  have  to  be  accounted  for 
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in  order  to  accurately  predict  the  turbulence  stress  component  via  a  seven  equation  Reynolds 
stress  transport  model.  The  rapid  response  to  the  local  flow  suggests  that  partially  resolved 
methods,  such  as  detached  eddy  simulations,  should  prove  useful,  where  only  the  shocks, 
expansions  and  flow  structures  shed  from  the  elements  need  be  resolved. 

The  strong  favorable  pressure  gradient  model  (Fig.  18b)  also  showed  large  periodic 
production  across  the  roughness  elements.  The  maximum  values  of  P  due  to  roughness  element 
leading  edge  compression  were  similar  to  those  observed  for  the  zero  pressure  gradient  models. 
However,  at  the  minima  within  the  trailing  edge  expansion  regions,  the  production  was  about  an 
order  of  magnitude  more  negative  that  for  the  zero  pressure  gradient  case  in  Fig.  17a.  Profiles  of 
the  allocation  of  {Px^i  and  {Pyy)x  had  the  appearance  of  complimentary  positive  and  negative 
cosecant  functions,  where  on  average,  75%  of  P  was  due  to  the  (Pxj:)/1  term.  Like  the  zero 
pressure  gradient  model,  in  the  leading  edge  compression  regions,  P  was  predominantly  made  up 
of  the  {Pxx)x  term,  while  it  was  equally  split  between  {Pxx)x  and  {Pyy)x  in  the  expansion.  The  {Pxy)x 
production  term  also  resembled  a  cosecant  function  that  was  in  phase  with  {Pyy)x-  On  average, 
{Pxy)x  was  about  -70%  of  the  value  of  P,  which  was  consistent  with  the  values  observed  in  the 
zero  pressure  gradient  smooth  wall  inner-region.  The  split  between  {Pxx)x  and  {Pyy)x  retained  the 
cosecant  appearance  exhibited  in  the  inner  region,  however  the  average  split  was  around  65/35 
for  both  of  the  outer  points  (y/S=  0.4  and  0.6).  The  distribution  of  (Pxy)x  was  similar  to  what  was 
observed  in  the  inner-region. 

5.3  Section  Summary 

The  near  wall  flow  structure  and  production  mechanisms  for  the  canonical  zero  pressure 
gradient  square  roughness  model  were  significantly  altered  from  the  conventional  bursting 
process  towards  localized  production  driven  by  roughness  element  induced  flow  structures  such 
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as  wake  vorticies  [Grass  (1971)].  The  net  effect  was  increased  turbulence  levels  across  the 
boundary  layer.  Near  the  wall,  the  principal  stress  axes  orientation  was  significantly  larger, 
which  resulted  in  a  redistribution  of  the  turbulence  stress  tensor  away  from  shear  and  towards 
compression.  Also  in  this  region,  the  production  tensor  was  highly  anisotropic,  where  the 
production  tensor  principal  axes  were  misaligned  with-respect-to  the  turbulence  stress  principal 
axes.  These  findings  explain  the  unequal  modification  of  the  turbulence  stresses  measured  in  a 
wind  axes  system. 

The  strong  favorable  pressure  gradient  axial  normal  stress  production  was  significantly 
reduced,  with  a  negative  value  neary/^=  0.4,  which  indicated  that  turbulence  energy  from  this 
component  was  being  transferred  back  to  the  mean  flow  as  suggested  by  Amette  et  al  (1996). 
Conversely,  the  transverse  normal  stress  production  was  positive  for  the  strong  pressure  gradient 
flow,  which  meant  that  this  component  was  directly  energized  by  the  mean  flow.  This  behavior 
is  opposite  to  that  of  the  zero  pressure  gradient  flow.  Thus,  the  overall  transport  dynamics  were 
significantly  altered  including  the  cross  coupling  between  the  axial  and  transverse  components 
(e.g.,  pressure-strain  redistribution).  Coupling  the  roughness  to  the  pressure  gradients  resulted  in 
a  competition  between  the  global  pressure  gradient  and  the  flow  structure  established  by  the 
roughness  topology. 

For  the  diamond  models,  the  production  was  periodic,  where  production  was  positive  and 
negative  on  the  forward  and  aft  half  of  the  roughness  elements,  respectively.  The  coupled  strong 
favorable  pressure  gradient  model  results  confirmed  that  shock  turbulence  interactions  were 
responsible  for  part  of  the  turbulence  level  increases  near  the  leading  edge  of  the  roughness 
element.  The  overall  flow  structure  and  turbulence  production  varied  significantly  on  a  local 
scale.  For  example,  the  principal  stress  axis  angle  exhibited  strong  variations,  between  -10°  to  - 
20°  in  the  favorable  and  adverse  pressure  gradient  regions,  respectively.  The  scalar  production 
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had  large  maximum  values  near  the  leading  edges  (similar  to  the  square  roughness),  where,  in 
the  expansion  regions,  negative  produetion  was  seen.  The  component  production  varied  across 
the  element,  where  the  axial  production  was  dominant  in  the  leading  edge  regions;  the  transverse 
was  dominant  in  the  areas  of  negative  production,  and  the  shear  stress  production  was  nearly 
constant  at  about  40%  of  scalar  production  across  the  element.  Coupling  the  diamond  and  strong 
favorable  pressure  gradient  had  the  net  effect  of  increasing  the  spatial  extent  of  the  favorable 
pressure  gradient  region,  which  resulted  in  larger  negative  values  of  the  scalar  production,  and  a 
redistribution  of  the  production  components. 
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(a)  Zero  Pressure  Gradient 

V/S^  0,4 


-JO 

1.00 

-20 

0.75 

-10-  ^ 

0.50 

0  ^ 

0.25 

10 

0.00 

h 

\ - ' - ’ - ' — ■ 

A 

F 

f*  ^  ^  A  ^ 

t  A  A  *  *  ^ 

A  ^  A  A  A  ' 

A  *  A  *  ‘ 

_ 

;  D  0  n  „ 

Q  n  n  n  □  1 

p""Lion[4]Lj- 

L. . 1 , , 

fcj  Lf  LI 

0*005 

0.004 

0.003 

0*002 

0*001 


0.000 


— . — 1 — . — . — 1 

'  ^  ■  1 

1 — ■ — ' — ■ — ' — : 

Id  ®  E 

a 

s  ®  E 

ta  -i 

^  E3  H  a  H 

'  Eja 

s  a  a  . 

29*0 


29*5 


30*0 

x/S 


30*5 


31*0 


ref 


y/S^0,6 


a 

r 


— 1  1  .  .  1 

1  .  .  .  . 

"  *  ^  i 

*  *  4 

^  *  * 

.  * 

A 

A. 

A  A 

*** 

‘  *  *  * 

]□□□□□□ 
— 1  .  . 1 

□  On  pocji; 

III.. 

1  □  □□• 

^  □  n^n 

0.004  . 
0.003 


0.00!  fHQg 


0.000 


-JO 

-20 

S 

-10- 


10 

0*005 


0.004 

0*003 


0.000 


15*0 


15*5 


16.0 

x/8 


16.5 


17*0 


ref 


y/S=0,6 


-JO 

1.00 

-20 

0.75 

h- 

-10-  ^ 

0.50 

0  * 

0.25 

10 

0.00 

0.005 

0*005 

0*0U4 

0*004 

0*003 

0.003 

0.002 

S  0*002 

0.001 

0*001 

0*000 

0*000 

■  1  .  .  . 

* 

1  .  .  ,  . 

* 

[-20 

*■  ^ 

J. 

A 

* 

-10-  ' 

r  *  i 

A  A  A  A  ^ 

[  ^  A  A  A  A 

p  o  a  □  □  n 

A  A  A  A  A  . 

n  n  n  n  n 

A  A  A  A  A  £ 

n  n  □  n  Q  li 

[ 

n  0  ' 

P  u  □  ii  n  ^ 

^ . J 

1 , , 

1  10 

0.005 

0*0U4 

1  0*003 

■  '  '  '  1 

1  ■  ■  '  ■  1 

1  ■  1  r  1 

i 

Sfl _ _ 

a  a  0  ^ 

.  El  ,0  ig  _  _ 

0.002 

0.001 

0*000 

13a®  . 

,*,  t  i,.j  ,j 

1  .  .  i  *  J 

29*0 


29*5 


30*0 

x/S 


30*5 


31*0 


ref 


(b)  Strong  Favorable  Pressure  Gradient 
Fig.  17  Axial  profiles  of  diamond  roughness  turbulence  stress 
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6.  Concluding  Remarks 


It  was  demonstrated  in  the  present  study  that  the  local  shock  and  expansion  waves  generated 
by  small  roughness  elements  (k/S~  0.05)  can  provide  a  mechanism  to  alter  the  mean  flow  field 
and  turbulence  structure  across  a  boundary  layer,  which  is  fundamentally  different  from  the 
subsonic  flow.  To  better  understand  the  response  of  the  mean  flow  and  turbulence  to  the 
roughness  elements,  a  series  of  experiments  were  performed  on  a  supersonic  (M  =  2.87)  high 
Reynolds  number  boundary  layer  (Ree  =  40,000).  The  two  roughness  topologies  consisted  of 
diamond  and  square  shaped  patterns.  To  provide  a  complete  and  relevant  basis  for  analyses, 
experiments  were  performed  under  the  same  conditions  and  with  the  same  instrumentation  for 
the  following  cases:  (1)  zero  pressure  gradient  boundary  layer  on  a  smooth  wall,  (2)  zero 
pressure  gradient  boundary  layer  with  square,  canonical  J-type,  roughness  topology,  (3)  zero 
pressure  gradient  boundary  layer  with  diamond  shaped  A:-type  roughness  topology,  which 
produces  a  strong  secondary  flow  structure  through  a  pattern  of  attached  oblique  shocks  and 
expansion  waves,  (4)  favorable  pressure  gradient  boundary  layers  on  a  smooth  wall,  (5) 
favorable  pressure  gradient  boundary  layers  with  square  roughness,  (6)  favorable  pressure 
gradient  boundary  layers  with  diamond  roughness,  (7)  combined  favorable  then  adverse  pressure 
gradient  boundary  layer  on  a  smooth  wall,  (8)  combined  pressure  gradient  boundary  layer  with 
square  roughness,  and  (9)  combined  pressure  gradient  boundary  layer  with  diamond  roughness. 

The  purpose  of  the  first  two  cases  was  to  provide  a  canonical  background  for  comparison 
purposes,  where  the  d-type  roughness  was  found  to  follow  the  subsonic  trends.  The  purpose  of 
case  three  was  to  characterize  interactions  between  the  local  pressure  waves  and  boundary  layer 
mean  and  turbulence  flow  structure.  Cases  four,  five  and  seven  through  nine  were  examined  to 
characterize  global  distortion  effects.  The  favorable  pressure  gradients  in  case  six  were  selected 
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to  isolate  and  demonstrate  the  importance  of  shock-turbulence  interactions  near  the  leading  edge 
of  the  roughness  elements.  The  heights  (0.8  mm)  and  the  frontal  blockage  for  both  topologies 
were  matched.  To  characterize  the  local  interactions,  the  experimental  data  were  acquired  with 
sufficient  resolution  to  characterize  the  mean  and  turbulence  flow  structure  over  individual 
roughness  elements.  Particle  image  velocimetry,  pressure  sensitive  paint,  pitot  pressure  profiles 
and  schlieren  photography  were  employed  to  achieve  this  experimental  goal.  An  extensive 
database  was  achieved  and  analyzed.  A  synopsis  of  the  resulting  basic  flow  processes  are  given 
below. 

The  diamond,  A:-type,  roughness  was  shown  to  alter  the  mean  flow  in  manner  that  was 
significantly  different  than  the  square,  <i-type,  roughness.  Specifically,  the  diamond  element 
shock  and  expansion  waves  had  a  measureable  effect  on  the  local  surface  pressure  and  mean 
flow  field.  The  square  roughness  topology  appeared  to  follow  the  sand  grain  behavior  described 
by  Goddard  (1959),  where  the  velocity  data  followed  the  incompressible  pattern  with  van  Driest 
II  (1951)  scaling.  When  spatially  averaged,  the  diamond  element  data  also  tended  to  collapse 
onto  the  law  of  the  wall  and  defect  law  plots.  These  results  indicated  that  the  average  behavior  of 
the  diamond  element  flow,  even  with  the  strong  wave  structure,  is  amenable  to  the  classical  low 
order  engineering  approximations.  The  zero  pressure  gradient  square  roughness  Reynolds 
stresses,  with  Morkovin  scaling  (1961),  for  the  smooth  and  square  models  collapsed  onto  the 
expected  trend  with  Townsend  scaling  (1976);  the  diamond  element  data  did  not.  The  magnitude 
of  the  square  roughness  Reynolds  shear  stresses  near  the  wall  was  nominally  2.4  times  that  for 
the  smooth  wall,  where,  for  the  diamond  model,  the  ratio  was  1.3.  This  implies  that  skin  friction 
was  significantly  lower  for  the  diamond  elements,  which  agrees  with  expectations  based  on 
Goddard’s  (1959)  model. 
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The  flow  stracture  and  production  mechanisms  for  the  canonical  zero  pressure  gradient 
square  roughness  models  were  shown  to  be  significantly  different  than  the  conventional  bursting 
process  typically  seen  for  smooth  wall  boundary  layers.  Near  the  wall,  the  principal  stress  axes 
orientation  was  found  to  be  significantly  larger,  which  resulted  in  a  redistribution  of  the 
turbulence  stress  tensor  away  from  shear  and  towards  compression.  Also  in  this  region,  the 
production  tensor  was  highly  anisotropic,  where  the  production  tensor  principal  axes  were 
misaligned  with-respect-to  the  turbulence  stress  principal  axes,  explaining  the  unequal 
modification  of  the  turbulence  stresses  measured  in  a  wind  axes  system.  Compressibility  effects 
were  found  to  be  minimal  and  the  flow  structure  across  the  boundary  layer  appeared  to  scale  with 
the  roughness  surface  condition  as  described  by  Townsend  (1976).  Hence,  prediction  of  the 
overall  trends  for  this  case  appears  to  be  amenable  to  current  two-equation  prediction  methods, 
where  the  roughness  is  treated  by  way  of  the  boundary  condition.  However,  the  redistribution  of 
the  production  tensor  indicated  a  modification  of  the  pressure-strain  redistribution,  which  is 
believed  to  be  important  for  accurate  prediction  of  the  turbulence  stress  components  via  a 
Reynolds  stress  transport  model;  this  implication  applies  to  all  of  the  distorted  canonical  flow 
cases. 

The  smooth  wall  strong  favorable  pressure  gradient  production  mechanism  was  significantly 
altered  compared  to  the  zero  pressure  gradient  flow.  Specifically,  the  strong  favorable  pressure 
gradient  axial  normal  stress  production  was  significantly  reduced,  with  a  negative  value  near  y/^ 
=  0.4,  which  indicated  that  turbulence  energy  from  this  component  was  being  transferred  back  to 
the  mean  flow.  Conversely,  the  transverse  normal  stress  production  was  positive  for  the  strong 
pressure  gradient  flow,  which  meant  that  this  component  was  directly  energized  by  the  mean 
flow.  This  behavior  is  opposite  to  that  of  the  zero  pressure  gradient  flow.  Thus,  the  overall 
transport  dynamics  were  significantly  altered,  including  the  cross  coupling  between  the  axial  and 
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transverse  components  (e.g.,  pressure-strain  redistribution).  The  Reynolds  stress  response  to  the 
local  strain  rates  followed  the  expected  trend,  where,  as  indicted  by  Bradshaw  (1974),  the 
distortions  affected  the  flow  by  an  order  of  magnitude  more  than  would  be  anticipated  by  the 
associated  extra  production.  This  result  suggests  that  current  two-equation,  or  higher,  turbulence 
models,  with  suitable  compressibility  corrections,  may  have  the  necessary  fidelity  to  capture  the 
correct  trends.  Coupling  the  square  roughness  to  the  pressure  gradient  flows  resulted  in  a 
competition  between  the  global  pressure  gradient  and  the  flow  structure  established  by  the 
roughness  topology.  The  competition  proceeded  in  an  expected  manner,  suggesting  that  this  flow 
may  also  amenable  to  current  prediction  techniques.  However,  in  addition  to  the  production 
tensor  redistribution  mentioned  in  the  previous  paragraph,  it  was  estimated  that  pressure-work 
terms  were  significant  in  the  turbulence  transport. 

The  diamond  roughness  element  was  found  to  produce  local  flow  phenomena  that  altered  the 
local  state  of  the  boundary  layer.  This  was  a  marked  departure  from  the  canonical  flows.  The 
overall  flow  structure  and  turbulence  production  was  found  to  vary  significantly  on  a  local  scale. 
For  example,  the  turbulent  stress  turning  angle  exhibited  strong  variations,  between  -10°  to  -20° 
in  the  favorable  and  adverse  pressure  gradient  regions,  respectively.  The  scalar  production  had 
large  maximum  values  near  the  leading  edges  (similar  to  the  square  roughness),  where,  in  the 
expansion  regions,  negative  production  was  seen.  The  component  production  varied  across  the 
element,  where  the  axial  production  was  dominant  in  the  leading  edge  regions;  the  transverse 
production  was  dominant  in  the  areas  of  negative  production,  and  the  shear  stress  production  was 
nearly  constant  at  about  40%  of  scalar  production  across  the  element.  Coupling  the  diamond  and 
strong  favorable  pressure  gradient  had  the  net  effect  of  increasing  the  spatial  extent  of  the 
favorable  pressure  gradient  region,  which  resulted  in  larger  negative  values  of  the  scalar 
production,  and  a  redistribution  of  the  production  components.  The  large  leading-edge 
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production  results  indicated  that  shock  turbulence  interactions  were  responsible  for  part  of  the 
turbulenee  level  increases  near  the  leading  edge  of  the  roughness  element. 

The  observed  intricate  Reynolds-stress/strain-rate  relationship  implied  that  the  turbulenee 
responded  strongly  to  the  loeal  distortions,  and  thus  numerical  simulations  must  resolve  the  flow 
structure  within  the  roughness  topology  to  accurately  capture  the  flow  structure.  The  rapid 
response  also  indicated  that  linear  effects  were  predominant  in  the  transport  dynamics.  The 
magnitude  of  the  local  Reynolds  stress  response  to  the  local  strain  rates  followed  the  patterns 
indieated  by  the  canonieal  flows,  whieh  suggests  that  eurrent  two-equation  turbulenee  models 
with  suitable  eompressibility  corrections  should  capture  the  local  trends  provided  the  flow  within 
the  elements  is  resolved.  The  observed  produetion  fields  indieated  an  involved  balanee  between 
production  and  redistribution,  and  this  balanee  responded  in  spatially  rapid  manner  to  the  local 
flow  distortions.  Again,  these  intrieate  processes  are  most  likely  important  for  accurate 
prediction  of  the  turbulence  stress  components  via  a  seven  equation  Reynolds  stress  transport 
model.  Lastly,  partially  resolved  methods  sueh  as  detached  eddy  simulation  may  be  suitable  for 
modeling  these  flows  with  strong  local  distortions,  since  these  methods  are  capable  of  resolving 
shocks,  expansions  and  flow  structures  shed  from  roughness  elements. 

This  work  was  sponsored  (in  part)  by  the  Air  Force  Office  of  Scientific  Research,  USAF, 
under  grant/contract  number  F49620-02- 1-0365.  The  views  and  conclusions  contained  herein 
are  those  of  the  authors  and  should  not  be  interpreted  as  necessarily  representing  the  official 
policies  or  endorsements,  either  expressed  or  implied,  of  the  Air  Force  Office  of  Scientific 
Research  or  the  U.S.  Government. 
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